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ARTICLE INFO ABSTRACT 

Article history: The subject of heat transfer enhancement has significant interest to develop the compact heat 

Received 23 September 2011 exchangers in order to obtain a high efficiency, low cost, light weight, and size as small as possible. 

Serre mae form Therefore, energy cost and environmental considerations are going on to encourage attempts to invent 
ay 


better performance over the existence designs. Streamwise vortices can be generated using small flow 
manipulators or protrusions such as wings and winglets configurations. Single-pair, single row, or two 
dimensional array of vortex generators (VGs) can be punched, mounted, attached or embedded 
Keywords: in the boundary layer of flow channel. VGs generate longitudinal and transverse vortices, while 
Heat transfer augmentation longitudinal vortices are more efficient for heat transfer enhancement than transverse vortices. A 
Vortex generators dramatic augmentation in thermal performance of the thermal system can be achieved but pressure 
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ieee š drop penalty is existed. Several parameters have been overviewed in this paper, which have 
pronounced effect on the convective heat transfer coefficient and pressure drop penalty. These 
parameters are: attack angle of VG, geometry of VG, standard and novel types of VG, spacing 
between the VG tips, number of pairs of VGs in the flow direction, rectangular or circular array 
arrangement of VGs, common-flow upper (CFU) or common-flow down (CFD) configuration of VG, 
pointing up (PU) or pointing down (PD) arrangement of VG with flow direction, Re number, channel 
aspect ratio, number of tubes of fin-tube heat exchanges (HE), circular or oval tubes of fin-tube HE, and 
location of VG respect to the tube of HE or from leading edge of the channel. This paper gives an 
overview about the early studies done in order to improve the performance of thermal systems with 
minimal pressure losses to derive systems with less negative impact on the environment and high level 
of energy economic. This study also provides an outlook for future work using nanofluids with vortex 
generators. 

This article is also summarizes the recent experimental and numerical developments on the thermal 
conductivity measurements of nanofluids, thermal conductivity enhancement, convection and con- 
duction heat transfer, some applications, main problems and suggestions for future works. 

© 2012 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Compact heat exchangers are widely used in such fields as 
automobile industries, power systems, heating and air condition- 
ing, chemical engineering, electronic chip cooling and aerospace, 
etc. This subject has received considerable attention to get high 
efficiency, low cost, small size, and light weight. For the vehicle 
industries, an additional benefit can be considered such as a small 
size for the radiator tends to get a small frontal area of vehicle 
which in turn leads to diminish the drag force and subsequently 
minimizes the fuel consumption. 

Several engineering applications had held an enormous atten- 
tion to improve heat transfer performance of the compact heat 
exchangers in gas-gas and gas-liquid side. For the gas, it is well 
known that the thermal conductivity of the gas is inherently 
lower than that of liquid needs to improve its convective heat 
transfer coefficient. The traditional heat exchangers are generally 
improved by serve enhancement techniques with concentration 
on many types of surface augmentation. The augmentation of 
heat transfer can be done by disrupting the boundary layer 
growth, increasing the turbulence intensity, and generating sec- 
ondary flows [1,2]. Fig. 1 shows several types of the compact heat 
exchangers. 

Longitudinal vortex generators (LVGs) are one of the passive 
methods. It can be mounted on channel surfaces to generate 
longitudinal vortices which create a secondary flow and disrupt 
the growth of the boundary layer, disturb the fluid flow and 
create fluid swirling. It in turn causes a heavy exchange of core 
and wall fluid and thereby improves the heat transport between 
the flowing fluid and channel walls. VGs have been extensively 
studied to their high heat transfer performance [3]. 

The enhancement of heat transfer by using VGs depends on 
many parameters such as their shape, their geometries, their 
position from the leading edge of the test duct, and angles of 


attack. The geometrical parameters of VGs have more effects on 
the heat transport. Four basic configurations of VGs are: delta 
wing, rectangular wing, delta—-winglet and rectangular winglet as 
shown in Fig. 2 [4]. The vortex generators VGs usually take 
configurations of tiny protrusions which may be incorporated 
into the main surface of the channel by punching, stamping, 
attaching, mounting and embossing [1]. Although the heat trans- 
fer surface area may be not changed before and after the set up of 
VG. The fluid flow can be significantly disturbed because of the 
longitudinal vortices generated by VG. Winglets have been 
successfully used in order to improve the convective heat transfer 
rate of modern thermal systems because they can generate 
intensive longitudinal vortices with less penalties in pressure 
drop [2]. 


2. Vortex generators 


The vortex may be divided into transverse vortex (TV) and 
longitudinal vortex (LV) according to its rotating axis direction. 
The axes of TVs lie perpendicular to the main flow direction, while 
LVs have their axes parallel to the main flow direction, thus they 
are also called streamwise vortices. In general, the LVs have 
been reported to be more effective than TVs on heat transfer 
enhancement [5]. A transverse vortices TV and longitudinal 
vortices LV are generated as a mainstream flows across the VG. 
The axis of the longitudinal vortices is parallel to the main flow 
direction whereas the axis of the transverse vortices is normal to 
the flow direction. Generally, the effect of longitudinal vortices 
is more effective than transverse vortices on heat transfer 


performance [1]. 
The fluid flow and heat transfer characteristics in flow chan- 


nels have been extensively studied in the past decades. This is due 
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Nomenclature 


C fanning friction coefficient 

Cp specific heat (J/kg K) 

d nanoparticle diameter (nm) 

f Darcy’s friction factor 

h heat transfer coefficient (W/m? K) 


k thermal conductivity (W/kg K) 
kg Boltzmann’s constant 

kej effective dielectric constant 
Nu Nusselt number, Nu=hD/k 
AP pressure drop, Pa 

Pe Peclet number, Pe=Lu/a 

Pr Prandtl number, Pr=C,u/k 
Re Reynolds number, Re=Du/v 
St Stanton number 

T temperature (K) 

u mean fluid velocity (m/s) 


Greek symbols 


attack angle 

,& volume fraction 
nanoparticle volume fraction 
density (kg/m?) 
dynamic viscosity (kg/m s) 
kinematic viscosity (m?/s) 


TEDOS DdD 


Subscript 

eff effective 

f base fluid 

fr freezing point of the base liquid 
in inner 

nf nanofluid 
np nanoparticle 
o base case 
out outer 

p particle 

s solid 


Abbreviation 


CFD common-flow down 

CFU common-flow up 

CHE compact heat exchanger 

CNT carbon nanotube 
C-MWNTcarbon multi-walled nanotubes 
ETC effective thermal conductivity 
LV longitudinal vortices 

PCM phase change material 

HTE heat transfer enhancement 
TV transverse vortices 

VG vortex generator 


to its wide applications in compact heat exchangers, solar 
collector, laser curtain seals, and cooling of electrical devices. 
However, in these applications, gases were used instead of liquids 
as heat transfer medium. Gases inherently have lower heat 
transfer rate than liquids. Therefore, heat transfer enhancement 
techniques that improve heat transfer rate are of particular 
interest. Many previous investigators have studied the techniques 
of heat transfer augmentation theoretically and experimentally. 
Heat transfer enhancement by VG is one passive method to 
generate streamwise vortices that create high turbulence in fluid 
flow over heat transfer surfaces. It was found that the only 
effective way to increase the heat transfer coefficient is to 
decrease the thermal resistance of the sublayer adjacent to the 


hot fluid 


cold fluid 


wall immediately where the viscous effects of the sublayer are 
dominant. This can be done by increasing the turbulence of the 
fluid flow in the main stream so that the turbulent eddies can 
penetrate deeper into this layer. The thermal resistance can be 
reduced and the heat transfer can be augmented by using VG that 
can be punched or attached to the internal surface of the channel 
wall. The vortex generator can be regarded as a special kind of 
extended surface, which can be stamped on or punched out from 
the fin. VG not only disturbs the flow field, disrupt the growth of 
the boundary layer, but also makes fluid swirling and causes a 
heavy exchange of core and wall fluid, leading to the enhance- 
ment of heat transfer between the flowing fluid and the channel 
walls. 


Fig. 1. (a) Plate-fin heat exchanger and its surface geometries, (b) plain rectangular fins, (c) plain triangular fins, (d) wavy fins, (e) offsets strip fins, (f) perforated fins and 


(g) louvered fins [3]. 
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—~ 


b/c 


A= 2 


Fig. 2. Longitudinal vortex generators types: (a) delta wing, (b) rectangular wing, 
(c) delta winglet pair and (d) rectangular winglet pair [30]. 


In the tube bank fin heat exchangers, the airside resistance 
generally comprises more than 90% of the total thermal resis- 
tance. Therefore, special surfaces are employed to effectively 
improve the airside heat transfer performance. VGs can be applied 
on the fin surface to generate longitudinal vortices in the main 
flow; hence the heat transfer of the fin surface is enhanced [6]. 

The techniques used to investigate the flow phenomena 
associated with VGs are classically IR imaging, particle image 
velocimetry (PIV), five hole pressure probes and the naphthalene 
sublimation technique which relates mass transfer to heat trans- 
fer by evoking the heat and mass transfer analogy. Hernon et al. 
[7] studied the phenomena associated with the VGs using the 
technique of hotwire and PIV images and made a comparison 
between them. They found some of the key observations from the 
instantaneous PIV are that at any given time the vortices shed 
from each delta winglet can vary in size, strength, and location 
from the wall. These variations cause the two counter rotating 
vortices to sweep across the lower wall while at the same time 
they move towards and away from the wall causing compression 
and extension of the underlying boundary layer. All of these 
complex unsteady flow phenomena may results in enhanced heat 
transfer. 

In the development of recent industrial world, the reduction in 
the heat exchanger volume and augmentation in their perfor- 
mance for heat transfer had a great interest of many researchers. 
In addition, there are always economical factors need to be 
considered to reduce the costs of heat transfer process. The need 
for high performance thermal systems due to the increase in the 
pressure drop associated with the augmentation in heat transfer 
becomes necessary. 

Henze and Wolfersdorf [8] examined the effects of longitudi- 
nal vortices induced by tetrahedral VG. They highlighted that the 
highest VG shows the highest heat transfer enhancement. The 
convection heat transfer coefficient increases as Re number 
increases while the heat transfer enhancement related to the 
smooth channel flow is decreasing. They found that the effect of 
turbulence is less remarked for regions which are affected by 
vortices. 


3. Applications of vortex generator 


Heat exchangers are widely used in many industrial areas such 
as chemical engineering, radiator of car, locomotive powered by 
internal combustion engine, geothermal, fossil, process plant, 
power system, air-conditioning, refrigeration, internally cooled 
gas turbine blades, air-cooled solar collectors, electronic chip 
cooling, wastewater aeration tanks, chemical mixers, nuclear 
reactors cores, aerospace, and other engineering applications. 
An efficient heat exchanger in such systems could results in the 
lesser consumption of the energy resource, which provides 
both economical and environmental benefits. Increasing demands 


are being placed on heat exchanger performance for reasons 
of compactness, economy in manufacturing and operating costs, 
energy conservation and even for ecological reasons. The 
importance of these issues continues to motivate the study of 
enhancement techniques. Universal types of heat exchangers 
incorporating plate-fins or fin-circular/flat tubes are designed 
for the above industrial processes and systems, and for these 
types the flow is channeled between plates [3,9-13]. VGs are used 
with the louvered fins to provide advantages over continuous fin 
designs because each louver initiates new boundary layer growth 
producing higher convective heat transfer than that of a contin- 
uous fin [14]. 

According to 2002 estimate by the US Department of Energy 
(EIA) heating ventilation, air conditioning, and refrigeration 
(HVACR) systems accounted for 55% of the energy used in 
residential buildings and 45% in the commercial buildings. This 
usage amounts to an annual expenditure of roughly 32 trillion 
dollars for energy consumed by HVACR systems alone. In the 
liquid-to-air and phase-change heat exchangers typical to HVACR 
systems, the air-side thermal resistance is the largest single 
contributor to the overall thermal resistance. For example, the 
air side can comprise 75% of the thermal resistance in an 
evaporator and 95% in condenser for typical refrigeration applica- 
tions. With rising energy costs and new legislation aimed at 
efficiency and environmental protection, heat exchanger perfor- 
mance will continue to be very important. Achieving even a 10% 
performance enhancement in this type of heat exchanger geo- 
metry could have profound implications on the HVACR systems 
where it is commonly used [15]. 

The need of high performance of thermal systems in many 
engineering applications has stimulated enormous interesting to 
find various methods to improve heat transfer rate in the system. 
The conventional heat exchangers are generally improved by 
various augmentation techniques which focus on many types of 
surface augmentation. There are two enhancement technologies 
of convective heat transfer for compact heat exchangers. One is to 
extend heat transfer surface area like a fin, the other is to increase 
heat transfer coefficients between solid surface and fluid such as 
using turbulators as vortex generators. In cooling channel and 
channel heat exchanger design, rib, fin, wing turbulators are often 
employed in order to increase the convective heat transfer rate 
leading to the compact heat exchanger and in turn leads to meet 
the desired of high efficiency and low cost with the volume as 
small as possible and the weight as light as possible. The use of 
wing turbulators completely results in the change of the flow field 
and hence the variation of the local convective heat transfer 
coefficient. Winglets have been successfully used for enhance- 
ment of heat transfer of modern thermal systems because they 
can generate longitudinal vortices and help to disrupt the main 
flow with less penalty of pressure loss [2,5,16-18]. 

Vortex induced heat transfer enhancement is a passive aug- 
mentation method that relies on the intentional generation of 
secondary, streamwise vortices to enhance heat transfer. The 
method may be utilized in many liquid-to-air and refrigerant- 
to-air heat exchanger applications. Streamwise vortices can be 
introduced using small flow manipulators such as those shown in 
Fig. 2. VG is used with the flat-plate compact heat exchangers to 
improve their thermal performance. 

For the array of VGs, regular arrays of longitudinal vortices 
embedded in turbulent boundary layer occur frequently behind 
evenly spaced pro-turbulences, and behind vortex generator rows 
installed to increase mixing. VG array embedded in turbulent 
boundary layer causes periodic variations in the skin friction with 
corresponding increase in heat transfer [19]. 

The most common practical application is the introduction of 
the vortex into a boundary layer to suppress separation. These 
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flows also occur naturally when horseshoe vortices form near the 
junction of a blunt body and a wall on which a turbulent 
boundary layer is developing. A typical example is the junction 
between a gas-turbine blade and its endwall, where the enhanced 
mixing due to the vortices can lead to locally high heat transfer 
rates along the endwall [20]. 

The subject of heat transfer is of serious interest in heat 
exchanger applications. The process industries call for more and 
more compact designs of heat exchangers. Flow interruption created 
in flow at periodic intervals is a popular means for heat transfer 
enhancement in compact heat exchangers. One relevant application 
using such flow configurations is the heat transfer between the 
flowing fluid and plates in the case of plate-fin heat exchangers [21]. 
Since global energy consumption is starting to impact negatively on 
the environment, as well as causing depletion of the existing fuel 
stocks, energy utilization has come under close scrutiny. Energy cost 
and environmental considerations continue to motivate attempts to 
derive better performance over the existing designs [22]. 


3.1. Circular and non-circular ducts 


Li et al. [23] illustrated that the vortex intensity and vortex 
flow length are strongly affected by the rib dimensions in laminar 
flow of water tube. 

Kenan et al. [24] studied experimentally the tapes with 
double-sided delta-winglets under different geometrical and flow 
parameters in a test pipe. They pointed out that the winglet 
height has the most effective parameter on the friction factor, 
while the most effective parameter affects the amplitude of the 
pressure fluctuations is the Re number. 

Several studies have been performed to study the influence of 
pair/row of VG mounted on a rectangular channel wall on heat 
transfer performance. Pauley and Eaton [19] studied array of one 
transverse row of VG with several spacing between the edge of 
VG which varied from 2 cm to 14 cm and different angle of attack 
which varied from 6° to 24°. They have found that one row of VG 
causes a tremendous heat transfer augmentation over large areas. 
They found that the peak augmentation in the Stanton number 
(St) was approximately 30% and the minimum value was 8% 
below the undisturbed level for all cases with different VG 
spacing. They observed that the downstream where the vortices 
have lifted themselves out of the boundary layer, the peaks in 
heat transfer augmentation became less pronounced and 7-18% 
level of heat transfer augmentation was seen across the channel. 
They pointed out that the co-rotating array produces 13% aug- 
mentation in heat transfer over the entire surface which is 
somewhat higher than observed for the alternating array. Wro- 
blewski and Eibeck [20] found that a peak enhancement in heat 
transfer performance was 25% in the region near the downwash 
side of DWVG mounted in a turbulent boundary layer when the 
attack angle is 12°. Kim and Yang [25] found that the common 


a 


-: 0.2 m/s 


Y (em) 


flow-down (CFD) cases showed better heat transfer characteris- 
tics than the common flow-up (CFU) cases when a pair of DWVG 
was embedded in turbulent boundary layer of rectangular chan- 
nel. They obtained, at f=45°, an enhancement in heat transfer 
rate for the CFD by 7% and 3% more than those of 6=20° and 30° 
respectively. The secondary velocity vectors at 8=30° for the CFD 
and CFU is shown in Fig. 3. 

Depaiwa et al. [2] examined the effect of existence of ten pairs 
of rectangular winglet vortex generators (RWVG) at the entrance 
of a rectangular solar air heater channel. The attack angle of VGs 
was varied from 30° to 60° pointing upstream (PU) and pointing 
downstream (PD) of the flow while Reynolds number was varied 
from 5000 to 23,000. They found that the larger attack angle leads 
to higher heat transfer rate and flow loss than the lower one. They 
displayed that the PD-RWVG yields higher heat transfer rate than 
PU-RWVG. They reported that convective heat transfer coefficient 
is strongly affected by the attack angle of VG. 

Qiuwang et al. [26] achieved experiments on the heat transfer 
and pressure drop in horizontal narrow rectangular channels with 
mounted LVGs. Two and three rows of VG pairs were examined 
with different arrangements. They found that the LVGs could 
greatly improve the heat transfer rate by 10-45% and the heat 
transfer performance of channel with LVGs on two sides are 
better than those on one side. 

Jian et al. [27] performed tests to investigate fluid flow and 
convective heat transfer characteristics in a narrow rectangular 
channel with and without four pairs of LVGs. For laminar flow 
regime, VG causes heat transfer enhancement of about 100.9% 
and flow resistance increase of only 11.4%., while for turbulent 
flow regime, VG causes heat transfer enhancement of above 87.1% 
and flow resistance increase of 100.3%. They concluded that LV 
can obviously enhance heat transfer of single-phase water, and 
increase flow resistance mildly. 

Min et al. [3] tested the influence of RWPVG attached on the 
rectangular channel wall on heat transfer characteristics. Re number 
was ranged from 5000 to 17,500 with various valued of f. They 
modified the shape of rectangular wing MRWP as shown in Fig. 4. 
Their results revealed that the average Nu number increases with the 


p R 
©) 
qi 


Fig. 4. Schematic view of modified VG: (a) layout of CFD wings and (b) modified 
rectangular wing [3]. 
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Y (em) 


Fig. 3. (a) Secondary velocity vectors with CFD and (b) secondary velocity vectors with CFU, at 6=30° [25]. 
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attack angle and 55° angle shows the highest average Nu number 
which was 46-55% compared with plain channel in the range of their 
study. MRWP have higher heat transfer enhancement and lower 
friction factors in comparison with RWP. 

Promvonge et al. [17] and Chompookham et al. [18] performed 
tests to investigate the heat transfer characteristics in a rectan- 
gular channel fitted with combined wedge ribs turbulators and 
winglet type VGs for the turbulent regime. Re number was varied 
from 5000 to 22,000. The results of Promvonge et al. [17] revealed 
that the combined ribs and DWVG give significant increase in 
heat transfer rate over the smooth channel. The increase in Nu 
number is around 240% over the smooth channel when in-line 
ribs and WVG are used. The Nu number obtained from combined 
ribs and DW pair of VGs was found to be much higher than those 
from the ribs only. This is around 50% greater than that of the rib 
alone. The rib along with the WVGs of higher attack angle value 
provides higher heat transfer than that of lower attack angle 
value. The larger the attack angle value leads to higher friction 
loss than the lower one. The friction factor of the combined WVG 
and ribs was found to be higher than that of the rib alone around 
25-125%. While the results of Chompookham et al. [18] found 
that the combined staggered rib and the WVGs with attack angle 
of 30° should be applied instead of using the rib/WVGs alone to 
obtain higher heat transfer performance of about 40-65%. The 
friction factor values obtained from combined the ribs and WVGs 
are found to be much higher than those from the rib or WVGs 
alone. 

Promvonge et al. [28] showed that the larger attack angle of 
the DW leads to higher heat transfer and friction loss than the 
lower one. They observed that the largest attack angle (60°) of the 
PU-DW gives the highest Nu number and friction factor while the 
lowest attack angle of the PD-DW yields the best thermal 
performance. 

Aris et al. [29,30] investigated experimentally and numerically 
the use of wing type surface protrusions in the form of 3-D delta 
wing tabs adhered to the fin surface, thin wings punched-out of 
the fin material and TiNi shape memory alloy delta wings which 
changed their angles of attack based on the fin surface tempera- 
ture. The 3-D tabs in a staggered arrangement gave an augmenta- 
tion in the heat transfer by about 37% compared to a plain fin 
stack. While the punched-out delta wings in the staggered and 
inline arrangements provided enhancements of 30% and 26%, 
respectively. The pressure drop through the fin stack was by 
approximately 19% in the de-activate position. Their CFD results 
were in good agreement with their experimental work and 
previous reports. While Aris et al. [30] obtained maximum heat 
transfer improvements of up to 90% and 80% by the single and 
double wings, respectively, in the laminar-transition air flow. The 
pressure drop penalty across the test section, when the wings 
were activated, increased between 7% and 63% of the losses at 
their de-activated positions, for the single and double VG, 
respectively. 

Mochizuki et al. [31] studied the fluid flow characteristics in a 
rectangular channel with DWPVG mounted in a turbulent flow 
regime. They reported that the mean velocity profile averaged 
over the spanwise direction has acceleration near the wall and 
deceleration away from the wall. 

Numerically, Charbel et al. [32] performed a 3-D simulation to 
study the heat transfer in turbulent vortical flows. Vorticity is 
generated by inclined vortex generators in a turbulent circular 
pipe flow with different configurations that fall into three cate- 
gories which are, rows of trapezoidal vortex generators in 
different arrangements; the vortex generators are fixed at certain 
distance from the tube wall, and vortex generator rows between 
which a row of small protrusions are inserted on the tube wall. 
They pointed out that the longitudinal variation of local Nu 


number computed on different high-efficiency vortex mixer cross 
sections shows the advantage of using protrusions between two 
successive vortex generator rows in which they increase the local 
heat transfer by increasing the temperature gradients and vorti- 
city very close to the heated wall. Saraç and Bali [33] performed 
experiments to investigate the heat transfer and pressure drop 
characteristics of a decaying swirl flow in a horizontal pipe. The 
decaying swirl flow is produced by inserting VGs with propeller- 
type for Re number ranged from 5000 to 30,000. Their results 
indicate that the Nu number increases ranged from 18.1% to 163% 
depending on Re number, the position of VG, the angle of attack 
and the number of vane. They observed a decrease in Nu number 
when the angle of attack increases. They attributed that to an 
increase in the vane angle would increase the swirl period. 

Wu and Tao [5,34] investigated the heat transfer character- 
istics in a rectangular channel with a pair of RWVG mounted on 
the lower wall of the channel. Re number was ranged from 800 to 
3000. They remarked that the thickness of RWVG can cause less 
heat transfer enhancement in the region near to the VG. They 
found that the pressure drop in the rectangular channel with 
RWVG increases rapidly with the increase of the attack angle of 
the LVG. They observed that the thickness of LVG has little 
influence on average friction factor of the channel. The LVGs can 
improve the synergy between the velocity and temperature fields 
not only in the region near LVG but also in the large downstream 
region of LVG. LVs enable to enhance the global heat transfer of 
the channel. Whereas TVs can only enhance the local heat transfer 
of the channel. The attack angle of 45° of VG always provides 
better effectiveness of heat transfer enhancement. They found 
that Nu number decreases with the LVG’s location away from the 
inlet of the channel and decreases too with the space between the 
LVG pair decreases. They have seen that Nu number increased 
with the increasing the area of LVG pair and the results of Nu 
number at S=45° is slightly higher than that of 6=60°. As the 
attack angle of DWVG increases from 15° to 60°, the convective 
heat transfer coefficient increases from 8-11% to 21-34%. They 
observed that the heat transfer was improved better with the 
increasing of the length of RWVG and decreasing the height of 
RWVG when the area of LVG was fixed. They found that the 
DWVG is more effective than RWVG on Nu number at the same 
area of LVG. The geometry of VG plays an important role in the 
enhancement of heat transfer [34]. 

Yang et al. [35] studied the effects of the common-flow-up pair 
produced by vortex generators in a rectangular channel flow on 
fluid flow and heat transfer. They displayed that the distortion of 
thermal boundary layer is not as strong as the distortion of 
hydraulic boundary layer. The influence of the common-flow-up 
pair on heat transfer enhancement is maintained at the down- 
stream location 30 times as large as the chord length of the vortex 
generators. 

Biswas and Chattopadhyay [36] solved numerically the com- 
plete Navier-Stokes and energy equations and have seen that the 
delta-wing VG has a strong effect on the heat transfer perfor- 
mance when punched on the rectangular channel wall. They 
observed that the heat transfer is enhanced by about 45.4% over 
the smooth channel case at a non-dimensional axial distance of 
4 from the inlet and at attack angle of 20°. They found that the 
improvement was 26.4% in Nu number for f/=30° over the case of 
f=20° at the same location. At the exit of the channel and 
Re=1815, Nu number was 98.38% greater than that when 
Re=500. Zhu et al. [37] studied the heat transfer performance 
for fully developed flow in a rectangular channel using RWVG on 
one wall and rib-roughness elements on the other wall. They used 
the SOLA-algorithm to solve the time-dependent Navier-Stokes 
equations directly for the primitive variables by advancing the 
solution explicitly in time. They displayed that Nu /Nuo was 3.42 
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when Re=10° and 8=45°. Maximum augmentation obtained by 
Zhu et al. [37] was 450% using this arrangement of VG and ribs. 
The f/f, is 32.0 at Re number equal to 15x 10* and B=45°. 
Ahmed et al. [12] examined several types of VG embedded in 
laminar and turbulent boundary layer of rectangular channel. 
They stated that the delta-wing shows the generation of counter 
rotating longitudinal vortices along their leading edges that 
diverge slightly depending on the aspect ratio of the wing. They 
have seen that the longitudinal vortices follow somewhat the 
wing inclination angle until they are redirected by the duct wall. 
Kaniewski et al. [38] found that one pair of RWVG attached on 
lower wall of rectangular channel can increase the global mixing 
by more than 50% for laminar flow. 

Riitten and Krenkel [39] tested the effect of existence of delta- 
winglet pair punched on the upper heated wall of rectangular 
channel in a flow form of CFU. While another pair of VG was 
punched on the lower cold wall in a flow form of CFD (facing one 
another). They stated that this arrangement of VG enhances the 
mixing of cold and warm fluid in the channel. The configuration of 
the VG in the channel and the temperature distribution is shown 
in Fig. 5. 

While Hiravennavar et al. [40] studied the influence of the 
existence of DWVG pair of non-zero thickness embedded in the 
hydrodynamically developed and thermally developing laminar 
rectangular channel flow. They have found that heat transfer 
enhancement increases by 33% when single winglet is used and 
by 67% when a winglet pair is utilized compared to the plain 
channel. It was observed that as Re number increases from 790 to 
1580 the heat transfer increases by 37.6% at the end of the 
channel. They obtained an increase in Nu number by about 
0.83-12.49% when W/H (thickness of VG/height of the channel) 
varied from 0.0622 to 0.2485. Biswas et al. [41] examined the 
presence of delta-wing and delta-winglet pair VGs on heat 
transfer performance embedded in laminar boundary layer of 
rectangular channel. This geometric configuration of the channel 
is considered from their viewpoint as a representative of single 
element of either a compact gas-liquid fin-tube crossflow heat 
exchange or a plate-fin crossflow heat exchanger. They have seen 
that the delta-wing VG is more effective than the winglet-pair but 
the use of winglets appears to be a more attractive augmentation 
technique. They reported that the channel may not be allowed to 
be so long (X > 8.0). They stated that another row of VGs may be 
mounted before this length. They found that the combined 
spanwise average friction coefficient due to the winglet pair is 
less than that due to the wing. 

Deb et al. [42] have made an effort to investigate the heat 
transfer characteristics in laminar and turbulent flows through a 
rectangular channel containing built-in DWVGs. This channel is 
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considered as a representative of single element of a plate-fin 
heat exchanger. They found that a significant enhancement is 
observed for the cases with built-in DWVGs over the correspond- 
ing cases without any obstacle. For example, they found that the 
enhancement is about 16% at the exit of the channel when Re 
number equal to 5000. They also obtained a turbulent intensity by 
about 10% at the inlet of the channel for the cases of Re=5000 and 
Re=15,000, 6=12° and with aspect ratio of 1.6. The effect of the 
Re number on the distribution of combined spanwise-average- 
skin-friction in the channel with built-in winglet pair is shown in 
Fig. 6. 

Sohankar [43] investigated the heat transfer augmentation ina 
rectangular channel with two angled ribs as a V-shaped VG. The 
VGs are attached on the bottom wall of the channel and their 
angles in respect to the main flow direction are between 10° and 
30°. Re number is ranged from 200 to 2000. They observed that 
the heat transfer becomes steady at lower Re number while it 
becomes unsteady at higher Re number values. They found that 
the vortex strength increases with Re number and it decays in the 
streamwise direction. The peak values of Nu number they 
observed occurs on the lower wall upstream of the VGs where 
the horseshoe vortices form. 

Munish et al. [44,45] stated that in order to obtain the heat 
transfer augmentation with RWVG the price of increase in 
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Fig. 6. Effect of the Re number on the distribution of combined spanwise-average- 
skin-friction in the channel with built-in winglet pair [42]. 
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Fig. 5. Vortices generated by VGs and their impact on the temperature field [39]. 
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pressure drop has to be paid. The pumping power increases with 
the increase of the attack angle of VG. A higher angle of attack a 
higher pressure drop. The pressure drop increases from 18.7% to 
31.3% when the attack angle of VG increases from 20° to 37° 
compared to the plain duct. For the constant area of VG, the 
pressure drop decreases from 23% to 18.7% with an increase in 
length from 0.894 to 1.23. 


3.2. Flat plate-fin heat exchanger 


A plate-fin heat exchanger with triangular fins is a good choice 
because of excellent compactness and cost effectiveness when 
compared to traditional shell and tube heat exchangers or plate 
heat exchangers. The performance can be further increased by 
using LVGs that produce longitudinal vortices that strongly 
enhance the heat transfer by mixing the fluid from the wall into 
the free stream and vice-versa. 

Nakod et al. [46] studied experimentally the effect of the 
finned surfaces and surfaces with VGs on the local heat transfer 
coefficient between impinging circular air jet and flat plate as 
shown in Fig. 7. VGs are in the form of equilateral triangle 
arranged in many circular rows. They observed an augmentation 
in heat transfer higher than that of the finned surfaces. This 
augmentation was as high as 110% for a single row of six VGs and 
Re=25,000 compared to the smooth surface. They reported that 
the effect of increase in the number of rows of VGs, arrangement 
of VGs in a row, radius of pitch circle of VGs and number of VGs in 
a row is to decrease the Nusselt number compared to the single 
row of six VGs. They found that the optimum configuration of the 
VGs is single row of six VGs. The inclination angle does not affect 
the heat transfer in case of the VGs. 

Joardar and Jacobi [47] investigated experimentally the effec- 
tiveness of 1500 delta-wings VGs mounted in full-scale wind- 
tunnel testing of a compact heat exchanger typical to those used 
in automotive systems. They reported that an average heat 
transfer increases over the baseline case of 21% for dry surface 
conditions and 23.4% for wet surface conditions with pressure 
drop penalty smaller than 7%. Gentry and Jacobi [48] studied the 
effect of delta-wing VG on heat transfer embedded in a laminar 
boundary layer of a flat-plate. Reynolds number was ranged from 
600 to 1000 and the attack angle of VG was from 10° to 55°. They 
found that the heat transfer was enhanced up to 50% compared to 
the plain flat-plate. They reported that the VG geometry and angle 
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of attack play important roles to predict how and where the 
vortex is started and how its position relative to the boundary 
layer evolves as it flows down the plate. 

Chen and Shu [49] showed that the presence of delta-wing 
vortex generator mounted on heated plate installed on the 
bottom wall of a duct in fan flows has little overall effect on the 
near-wall averaged axial mean velocity and axial vorticity, but 
increases the turbulent kinetic energy. They concluded that the 
delta-wing vortex generator in fan flows has little effect on the 
heat transfer augmentation. 

Kotcioglu et al. [21] and Tiggelbeck et al. [50] investigated the 
heat transfer characteristics in a plate fins channel using the four 
basic types of VG and RWP of VG respectively. Kotcioglu and his 
group studied the arrangement of the wings as an array of rows 
with attack angle which was varied from 7° to 20° positively and 
negatively with the main flow direction in a rectangular channel. 
They pointed out that the mixing effect in the intermediate region 
between wing cascades improves the heat transfer characteristics 
as the inclination angle of winglet increases. They indicated that 
the Nu number increases with the Re number and the dependence 
of Re number on the friction factor is strong. The results of 
Tiggelbeck et al. [50] imply that the winglets give better perfor- 
mance in heat transfer than the wings and a DWP performs 
slightly better than a RWP at higher attack angles greater than 
30°. They showed that DWPVG enhanced heat transfer by 46% at 
Re=2000 and 120% at Re=8000 over the heat transfer on a plate. 
Their results showed that the maximum heat transfer for the case 
with RWP of VG occurs at a f between 45° and 65° and this 
maximum is smaller than that of DWP of VG. While at f=90° a 
sharp decrease in Nu number was observed. Nu number increases 
almost linearly with Re number as shown in Fig. 8. The friction 
coefficient for channel with VGs starts slowly increasing whereas 
the friction coefficient for a smooth channel still decreases. 

Wu and Tao [51] made an experimental and numerical study 
to investigate the average convective heat transfer on the top and 
bottom surfaces of a plain plate and four plates with a pair of 
delta winglet longitudinal vortex generator punched directly from 
the plates with attack angles ranged from 15° to 60°. The average 
Nusselt number of the plate with attack angle of 60° is slightly 
higher than that of plate with attack angle of 45° yet may bring 
larger pressure drop. They obtained a deviation of 10% between 
their numerical data and experimental results which validates 
their numerical model. 
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Fig. 7. Details of vortex generators [46]. (a) Arrangement of single row of six vortex generators. (b) Vortex generator. (c) Arrangement of two rows of vortex generators 
(same circumferential pitch). (d) Two rows of vortex generators (inline arrangement). (e) Two rows of vortex generators (staggered arrangement). (f) Arrangement of single 


row of three vortex generators. 
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Fig. 8. (a) Complete fin area-averaged Nu number, Nu, normalized with Nu, and 
(b) the friction coefficient normalized with Cro as function of £ [51]. 


Experiments to study the enhancement of heat transfer in 
louvered fin heat exchanger using winglets placed on louvered was 
performed by Sanders and Thole [14]. They tested the effect of attack 
angle, aspect ratio, direction, and the shape of winglet in a laminar 
regime. They highlighted that the best heat transfer augmentation is 
found with RW pair of VG and it is around 3%, 36% and 38% at Re 
numbers of 230, 615 and 1016, respectively. They reported that in an 
attempt to optimize these winglet parameters, heat transfer aug- 
mentations as high as 39% were achieved with associated friction 
factor increment as high as 23%. The maximum augmentation in heat 
transfer obtained is 25% when the winglet direction and orientation 
were alternated on every other louver, while the maximum value is 
around 33% when all winglets were placed in the VG-B orientation 
with alternating direction as shown in Fig. 9. The RW pair of VG 
yields higher friction factor values than the VG-B winglet. The friction 
factor significantly increases as the angle of attack and the winglet 
size increase. 

Guobing and Qiuling [52] performed a new shape of VG of 
curved trapezoidal winglet CTW and compared with traditional 
vortex generators; rectangular winglet, trapezoidal winglet and 
delta winglet. They showed that DWP is the best in laminar and 
transitional flow region, while CTWP has the best thermohydrau- 
lic performance in fully turbulent region due to the streamlined 
configuration and then the low pressure drop, which indicates the 
advantages of using this kind of VGs for heat transfer enhance- 
ment. The CTWP showed that smaller attack angle, larger curva- 
ture and larger angle of inclination gives better thermohydraulic 
performance under the same conditions. An appropriate spacing 
between the leading edges of CTWP should be considered for 
different flow regions. They concluded that double rows of CTWP 
do not show better thermohydraulic performance due to the 
larger pressure drop and the spacing between the two rows of 
CTWP should also be optimized. 
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Fig. 9. Side, top, and bottom view of VG winglet [14]. 


Numerically, Li-Ting et al. [53] reported that the flat-plate 
channel with DWP has better overall performance than RWP. The 
CFD and CFU configurations of DWP have almost the same overall 
performance. The CFD configuration has a better overall perfor- 
mance than the CFU configuration for RWP. 

Munish et al. [44,45] and Akcayoglu [54] performed their 
investigation about the HTE in a plate-fin heat exchanger with 
triangular fins as insert using a winglet pair of VG. They observed 
an enhancement in heat transfer when a RW type VG is mounted 
on triangular fins. The attack angle of VG was varied from 20° to 
37° and Re number was ranged from 100 to 200. A constant wall 
temperature is considered as a thermal boundary condition and 
the air is taken as a working fluid. They confirmed that the 
average Nu number increases with an increase in the angle of 
attack and Re number. For the same area of LVGs, the increase in 
the length of the LVGs brings more HTE than increasing the height 
[44]. While in [45], the HTE was 13% over the case of without 
winglet pair even at the exit when Re=200, f=20°. They 
observed that the performance of the winglet pair increases as 
Re number increases from 200 to 500 and increases with the 
increase in the height of the winglet pair. Akcayoglu [54] studied 
the effect of double rows DW pair of VG in CFD and CFU 
configuration mounted on the duct’s slant surface. They found 
that the second flow is stronger behind the second pair of VG than 
behind the first pair but becomes weaker far from the second pair 
in the case of CFD. They observed that both configurations are 
able to create a counter-rotation. They showed a large vortex 
formation area and a great induced vorticity field between vortex 
pairs for CFU arrangement compared to the CFD arrangement. 

Ferrouillat et al. [55] studied the effect of DWP and RWP VG on 
the performance of compact heat exchanger. They stated that for 
different turbulent models used, that DWP is more efficient than 
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RWP in terms of compactness criterion. In addition, they found 
that the optimum heat transfer and mixing efficiencies can be 
obtained when the distance between VG rows is around 7-10 
times the channel heights. 

Ahmed and Lars [56] showed that the temperature fluctua- 
tions, turbulent kinetic energy and unsteadiness effects are 
stronger in the region of a plate-fin heat exchanger where the 
longitudinal vortices are more active. 

Leu et al. [57] performed a combined numerical and experi- 
mental study of the heat transfer and fluid flow characteristics in 
the 3-row plate-fin and circular tube heat exchanger with 
inclined block shape VG mounted behind the tubes. They stated 
that the proposed heat transfer enhancement technique is able to 
generate longitudinal vortices and to improve heat transfer 
performance in the wake region. They pointed out that this kind 
of VGs not only produces longitudinal vortices, but also aid the 
fluid into the wake re-circulation zone. They found that as the 
attack angle increased, the strength of the longitudinal vortex 
intensified and both the Colburn and friction factors increased. 
They demonstrated that the case of f=45° provides the best heat 
transfer augmentation with 8-30% increase in Colburn factor 
across the range of Re from 400 to 3000, while the Fanning 
friction factor is only increased by 11-15%. This gave the greatest 
area reduction ratio by 25% as shown in Fig. 10. 


3.3. Tube-fin heat exchanger 


An augmentation of heat transfer in a compact plain-fin-and- 
tube heat exchanger with presence of winglet VGs have been 
studied by [4,9-11,16,22,58,59]. Torii et al. [9] found that the heat 
transfer was enhanced by 30% and 10% for Re number ranging 
from 350 to 2100 with the using of DWVGs of the common flow 
up CFU for staggered tube banks heat exchanger. While the 
augmentation was 20% and 10% for in-line tube banks for the 
same range of Re number. Zhang et al. [10] found that VGs should 
be mounted as near as possible to the tube wall in the case of 
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three-row flat tube bank fin, and vortices generated by the 
upstream VGs converge at wake region of flat tube. They have 
seen that the friction factor decreases as Re number increases. Yoo 
et al. [16] stated that the overall average heat transfer coefficient 
of fin-flat tube heat exchanger with DWVGs is increased by 75% 
compared to that of fin-flat tube without DWVGs and is increased 
by 45% compared to that of the fin-circular tube heat exchanger 
without DWVGs. They observed that the pressure drop of the fin- 
flat tube bank with VGs is increased by 80% compared to that of 
the fin-flat tube bank without VGs, but is increased by 50% 
compared to that of the fin-circular tube bank without VGs. 
Zhang et al. [10] studied the effect of the pitch of in-line delta- 
winglet VGs on local and average Nu number of the finned three- 
row flat tube bank. The Re number ranged from 500 to 5000. They 
reported the following: the pitch of in-line VG has a large effect 
on heat transfer, VGs with small pitch can enhance the heat 
transfer of the fin surface downstream, and the effect of the pitch 
on heat transfer performance depends on the attack angle of VG 
for identical mass flow rate. 

Allison and Dally [22] demonstrated that, when the VGs placed 
directly upstream of the tube of the fin and tube radiator, the 
winglet surface has 87% of heat transfer capacity of the louvered 
fin surface. They found that 39° attack angle of VG gives the best 
flow structure. While Jalil et al. [58] examined four different 
shapes of winglet VGs around heated cylinder with attack angle 
varied from 20° to 32° and Re number ranged from 7200 to 
14,400. They concluded that the heat transfer was augmented by 
about 14% when winglets were used and this enhancement 
increases with the Re number and the attack angle of VG while 
the pressure drop was slightly increased; the local pressure 
coefficient increases when the winglet is used; the local pressure 
drop increases as the attack angle of VG increases; the pressure 
drop is relatively proportional with the Re number. They reported 
that the position of VG plays an important role on the perfor- 
mance of heat transfer. They showed that the trapezoidal winglet 
gives a better enhancement in heat transfer rate compared to the 


Fig. 10. (a) Local Nu number distribution across the fin surface at Rep; = 3000 ((a) VGo, (b) VG30, (c) VGas), (b) flow visualization for three different span angles at Rep, =500 
((a) VG30, (b) VGas, (c) VGeo), and (c) flow visualization for span angle of 45° at three different Repr ((a) Repn=500, (b) Repn=1500, (c) Repn=2500) [57]. 
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Fig. 11. VG implementation: (a) winglets are manufactured as strips with six winglets, (b) VG in fin-spaces of the heat exchanger provides approximately 3500 winglet 
VGs in CFU configuration and (c) heat exchanger geometry showing actual evaporator [59]. 


rectangular, triangular and circular winglets. Joardar and Jacobi 
[59] examined delta-winglet pairs arranged as a single-row VG 
placed at the leading tube and three-row-inline VG array placed 
at alternate tubes as shown in Fig. 11. They found that the 3VG 
array yields higher pressure drop penalty which decreases rapidly 
with increase in Re number. They reported that when the heat 
transfer is augmented by using VG, the size of heat exchanger can 
be minimized to get the same heat duty. So the available size 
reduction could be redeemed by reducing the flow length which 
leads to reduce the pressure penalty and can reduce cost and 
allow more compactness. 

Kwak et al. [60] and Tang et al. [61] found that the increase of 
VG attack angle or length, or decrease of VG height may enhance 
the performance of vortex-generator fin of the fin-and-tube heat 
exchanger. 

Shi et al. [4] studied the optimal fin spacing for three-row flat 
tube bank fin mounted VGs. They observed that for commonly 
used fin materials and fin thickness, the optimal fin spacing is 
about 2 mm in industrial application. Wang et al. [11] showed that, 
for the same height of winglet, the delta-winglet exhibits intensely 
vertical motion and flow unsteadiness. They interpreted that this 
leads to a better flow mixing phenomenon. They showed that the 
corresponding pressure drops of the DW pair of VG are lower than 
those of the annular winglet at the same Re number. The pressure 
drop penalty of VG is relatively proportional with Re number. 

Numerically, many researchers have studied the heat transfer 
augmentation in inline/staggered circular tube heat exchangers 
with several types of VG mounted on the channel walls such as 
[13,16,42-48]. Gorji et al. [13] studied the heat transfer on 
staggered circular tube heat exchanger with VG embedded in 
hydrodynamically developed and thermal developing laminar 
boundary layer. For an incompressible fluid, they have seen that 
the average Nu number increases with presence of airfoil obstacle 
by about 20-30% and these obstacles are effective to increase the 
heat transfer. They reported that the more Re number the better 
effect on Nu number. Their results revealed that a triangular 
obstacle with three attack angles improve the heat transfer 
compared to the plain case. The downstream reattachment 
lengths were found to depend only on the triangle’s angle of 
attack [62]. Joardar and Jacobi [15] studied the effect of three 
arrangements of DWVG on heat transfer performance of seven- 
row inline-tube heat exchanger. These arrangements are: single 
pair, three inline array, and three staggered array mounted in 
laminar boundary layer. It was observed that the CFU arrange- 
ment of DWVG has high effectiveness on the heat transfer 
enhancement. At a constant tube-wall temperature, the 3-VG- 
inline-array configuration achieved enhancement up to 32% over 
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Fig. 12. Local heat flux distribution (z=0) at different axial locations on the VG fin 
for four configurations [15]. 


the baseline case. The local heat flux on the VG mounted fin (z=0) 
for four arrangements are shown in Fig. 12. 

Lei et al. [63] found, that the heat transfer coefficient increases 
as the attack angle of VG increases when Re number was ranged 
from 600 to 2600 and the attack angle was varied from 10° to 50°. 
The DWVG with an attack angle of 20° and aspect ratio of 
2 provides the best performance over the range of Re number. 
Their results implied that the VG placed in CFU orientation are 
more effective flow at lower Re numbers than at high Re numbers. 
They found that the heat transfer coefficient increases with the 
increase in Re number. The results show that the DWVG in CFU 
configuration can enhance heat transfer performance without an 
excessive amount of pressure drop penalty. This pressure drop 
penalty increases with the increase of attack angle and the 
increasing tendency is strengthened gradually. 

Gorji and Soleimani [64] revealed that heat exchangers with 
triangle obstacle improve the heat transfer in comparison with heat 
exchangers without obstacle. The upstream location of flow separa- 
tion and the downstream reattachment lengths are found to depend 
only on the triangle’s angle of attack. Ünal and Gören [65] performed 
3-D computational simulations of the flow around a circular cylinder 
fitted with VGs. Their numerical calculations emphasized the effec- 
tiveness and the performance enhancing character of the VGs. Li et al. 
[66] performed a 3-D numerical simulation on laminar heat transfer 
and flow characteristics of a slit fin-and-tube heat exchanger with 
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Fig. 13. (a) Schematic diagram of core region of a slit fin-and-tube heat exchanger with LVGs, (b) streamlines generated by LVGs and (c) temperature distribution in four 


y-z sections at Re=1051.3 [66]. 


LVGs. The core of the compact heat exchanger, the streamlines and 
the temperature distribution at four y-z sections with the using of 
LVGs is shown in Fig. 13. 

They found a significant heat transfer enhancement with the 
modest pressure drop penalty. The same study was performed by 
Zeng et al. [67] to examine the effect of DW pairs of VG on the 
thermal performance of fin-and-tube surface. They observed that 
the intensity of heat transfer can be greatly increased with the 
increase of the VG attack angle and the length and height of VG 
accompanying with the increase in pressure drop. Wu and Tao 
[68] displayed that the heat transfer enhancement of delta- 
winglet with an attack angle of 45° is larger than that of 30° 
while the pressure drop associated is also greater for the range of 
Re number of 800-2000. 

Tiwari et al. [69] and Chu et al. [70] found that the average Nu 
number for the three-row fin-and-oval-tube heat exchanger with 
LVGs increases by 13.6-32.9% over the baseline case, and the 
corresponding pressure loss increased by 29.2-40.6%. They dis- 
played that LVGs with placement of downstream, angles of attack 
30° and two tube-row provide the best heat transfer performance. 

Other researchers have investigated the heat transfer enhance- 
ment (HTE) in many-row oval (flat) tube heat exchangers using 
different types of VGs such as [6,71-73]. Chang et al. [6] 
examined the influence of the attack angle and the aspect ratio 
of a DWVG which is punched near the leading edge of the fin of a 
compact heat exchanger (CHE) on heat transfer enhancement. The 
VG was embedded in a thermally and hydrodynamically devel- 
oping laminar boundary layer with angle of attack ranged from 
20° to 45°. They displayed that the winglet with f=30° and 
aspect ratio of 2 provides the best ratio of HTE. They highlighted 
that when the dimensionless distance of DWP from the leading 
edge changes from 1.46 to 0.63, a little HTE has. Chen et al. [71] 
punched LVGs in form of winglets in staggered arrangements to 
enhance heat transfer in finned oval tube heat exchanger ele- 
ments at Re=300 and f=30°. They demonstrated that the 


staggered arrangement of the winglets has more effective than 
the inline arrangement for heat transfer enhancement. In the 
staggered arrangement, which is away from the tube, has more 
effective for heat transfer than winglets which is near the tube. 
They also found that the temperature distribution in the flow 
passage is intensively distorted by the staggered longitudinal 
vortices. The vector-plots of the secondary flow and correspond- 
ing 2D streamlines in eight cross-sections of several configura- 
tions of VGs at Re=300 are shown in Fig. 14. 

Chen et al. [73] added a new idea to augment the heat transfer 
in a flat-tube bank fin surfaces where VGs are mounted on both 
surfaces of the fin. Their results revealed that using this configura- 
tion of VGs, the height of VGs can be reduced and still obtain 
satisfactory THE. They found that if VGs on one surface of the fin is 
determined, the locations where VGs are mounted on other surface 
of the same fin are very important. While Song et al. [72] studied 
the relationship between the intensity of the secondary flow 
produced by VGs which are mounted on a three-row flat tube 
bank fin surfaces and the strength of convective heat transfer. The 
results revealed that the averaged absolute vorticity flux normal to 
the main flow direction can reflect the intensity of the secondary 
flow. In most cases, the secondary flow cannot change greatly the 
boundary layer characteristics especially for the internal flow at 
beginning region of boundary layer. The averaged absolute vorti- 
city flux can account for only the secondary flow effects on 
convective heat transfer but cannot quantify the effects of devel- 
oping boundary layer on convective heat transfer. 

Lemouedda et al. [74] performed a numerical investigation to 
study the performance of the finned-tubes using DWVGs. They 
demonstrated that VGs increase the level of vorticity especially 
behind the tubes and increase the mixing inside the finned tubes. 
This process enhances the heat transfer rate. They found that the 
tubes are exposed to stronger mean flow where the first pair of 
VG is in CFU configuration with mainstream flow compared with 
CFD configuration. This leads to increase the local heat transfer 
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Fig. 14. Secondary flow and the corresponding 2D streamlines on eight cross-sections in several configurations at Re=300: (a) velocity vectors and (b) temperature 
distribution [70]. 
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Fig. 15. Streamlines and temperature distributions at a horizontal section Z=0.065H at Re=600: (a) inline arrangement ((a) Without winglets, inline, instantaneous, 
(b) without winglets, inline, time-averaged, (c) angle of attack f=60°(common-flow-up), inline, time-averaged, (d) angle of attack p=+60°(common-flow-down), inline, 
time-averaged) and (b) staggered arrangement ((a) Without winglets, staggered, instantaneous, (b) without winglets, staggered, time-averaged, (c) angle of attack f=60° 
(common-flow-up), staggered, time-averaged, (d) angle of attack p=+60° (common-flow-down), staggered, time-averaged) [74]. 


accompanied by increase in pressure drop. They have seen that 
increase the number of DWVGs will increase the heat transfer 
rate which is also accompanied by increasing the pressure drop. 
The optimal design they obtained is 14% higher than that without 
VGs. The streamlines and temperature distributions at Z=0.065H 


and Re=600 can be seen in Fig. 15. They stated that this increases 
the flow resistance leading to higher pressure losses and hence to 
higher power input required. 

O'Brien et al. [75] showed that the addition of single pair of 
winglet to the oval-tube geometry yields significant heat transfer 
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enhancement averaged 38% higher than the oval-tube without VG. 
While the corresponding increase in friction factor was very modest, 
less than 10% at Rep,=500 and less than 5% at Rep»=5000. They 
demonstrated that the highest mean heat transfer coefficient can be 
observed for the case of circular-tube with winglets located on the 
downstream side of the cylinder, oriented at 45° angle to the flow. 


3.4. Heat sink 


Dake and Majdalani [76] designed thin winglet-type VG 
mounted on the leading edge of heat sink fins as shown in 
Fig. 16. This idea is applied for PCs, routers and servers. An 
experimental simulation of the interfin channel reveals that 
certain VGs, when placed upstream, can outperform others in 
their ability to fill the channel with pairs of strong vortices. 
Multiple pairs can also be generated to further accentuate the 
heat transfer using dual VGs. They summarized their results as 
follows: the flat shape of VG provides the most stable structure; 
the trapezoidal design induces the most oval shaped vortices with 
the spacing being a function of the width of the tip; the curved 
shape tends toward smaller; faster rotating vortices with the 
separating rotational center; the double delta design provides the 
least productive vortices exhibiting small outer diameter and 
widely spaced centers. Their results show that the lower angle of 
attack produces the widest vortex spread, presence of the com- 
mon-downwash region, and recirculation zones in each case. 

Kai-Shing et al. [77,78] stated that the enhancement using VG 
mounted in a heat sink are relatively effective when the flow is in 
the developing region whereas they become quite less effective in 
the fully developed region especially when the fin pitch is small 
or operated at a lower frontal velocity. They suggested that the 
vortex generators operated at a higher frontal velocity and at a 
larger fin pitch are more beneficial than that of plain fin geometry. 
The semi-circular VG has the highest heat transfer coefficients 
and pressure drops at developing region. They obtained a reduc- 
tion of 31.1% in the heat dissipation area at a frontal velocity of 
5 m/s within a developing region while for fully developed region 
the area reduction was of 1.8-11.5% at a frontal velocity 3-5 m/s. 
The results suggest that the vortex generators operated at a 
higher frontal velocity is more beneficial than that of plain fin 
geometry. The triangular attack vortex generator is regarded as 
the optimum enhancement design for it could reduce 12-15% 
surface area at a frontal velocity around 3-5 m/s. 


3.5. Electronic chips 


Investigations have been done to study the augmentation in 
heat transfer in air cooling of a staggered array of rectangular 
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Fig. 16. Shapes of vortex generators [76]. (a) Flat triangular design. ( 
dal design. (c) Curved design. (d) Double delta design. 
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electronics modules using DWVGs. Chomdee and Kiatsiriroat 
[79,80] observed that the VGs could enhance the heat transfer 
coefficient by 10-30% and reduce the thermal wake function and 
the module temperature when VGs integrated in front of all 
modules. They encouraged that the attack angle should not be 
more than 10°. They monitored that when VGs integrated in front 
of first row only, the heat transfer enhancement is distinct only at 
the first row. They observed that a higher attack angle results 
higher pressure drop in the fluid flow. They advised that the 
attack angle should not be more than 10°. 


3.6. Mini-channels and microchannels 


Qiuwang et al. [26] carried out experiments on the heat transfer 
and pressure drop in horizontal narrow rectangular channels with 
mounted LVGs for water flow to enhance heat transfer in cooling 
channels of plate-type fuel elements in reactor cores. They displayed 
that the LVGs could greatly improve the heat transfer rate by 
10-45%, and the heat transfer performance of channel with LVGs 
on two sides are better than those on one side. 

Jian et al. [27] stated that single-phase heat transfer of water can 
be enhanced greatly with acceptable increasing flow resistance after 
using several pairs of LVG as disturbing elements in narrow rectan- 
gular channel. For laminar flow regime, heat transfer augmentation 
was around 100.9% and flow resistance was increased only 11.4%, 
whereas heat transfer augmentation was over 87.1% and the flow 
resistance was increased 100.3% for the turbulent flow regime. They 
found the performance of heat transfer enhancement became better 
with increasing Re numbers in turbulent regime. 

Chao et al. [81] showed that heat transfer is enhanced in 
rectangular microchannel using LVGs compared with the smooth 
microchannel, with large pressure drop. They found that changing 
the direction of attack angle of one pair of LVGs in the streamwise 
direction can delay the occurrence of laminar-to-turbulent transi- 
tion to some extent. 


3.7. Refrigeration 


In most domestic and commercial refrigeration systems, frost 
forms on the air-side surface of the air-to-refrigerant heat exchanger. 
Frost-tolerant design typically employ a large fin spacing in order to 
delay the need for a defrost cycle. Unfortunately, this approach does 
not allow for a very high air-side heat transfer coefficient, and the 
performance of these heat exchangers is often air-side limited. 
Sommers and Jacobi [82] stated that LVG is an effective technique 
for thinning the thermal boundary layer and enhancing heat transfer, 
but its efficacy in a frosting environment is essentially unknown. 
They studied an array of delta-wing VGs applied to a plain-fin-and- 
tube heat exchanger as shown in Fig. 17. They indicated that the air- 
side thermal resistance is reduced by 35-42% when VG is used. 
While the heat transfer coefficient increases by 60-93% when Re 
number is ranged from 500 to 1300. 

The heat transfer performance of the air-cooled condensers used 
in binary geothermal power plants can be developed by replacing the 
circular tubes with oval tubes and adding strategically located VGs 
(winglets) in the fins [83]. They indicated that the addition of 
winglets increases the heat transfer coefficient by about 35% 
approximately, compared to plain tubes. Corresponding increase in 
friction factor is in the range of 5-10% for Re number based on the 
hydraulic diameter in the range of 500-5000. 


4. Nanofluids 


Advances in high power technology are demanding size reduction 
of thermal systems resulting in constrained spaces for increased heat 
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Fig. 17. (a) Heat exchanger geometry and VG configuration and (b) frost deposition on delta-wing type VGS [82]. 


transfer. This has led to thermal management problems in emerging 
systems which have defied solution by conventional cooling methods. 
Due to recent developments in nanotechnology, a new class of heat 
transfer fluids called nanofluids was discovered [84]. The nanofluid is 
a solid-liquid mixture in which metallic, nonmetallic or polymeric 
nanoparticles are suspended. The suspended ultrafine particles 
change transport properties and heat transfer performance of the 
nanofluid, which exhibits a great potential in enhancing heat transfer. 
However, the usual slurries, with suspended particles in the order of 
millimeters or even micrometers may cause some severe problems. 
The abrasive action of the particles causes the clogging of flow 
channels, erosion of pipelines and their momentum transfers into 
an increase in pressure drop in practical applications. Furthermore, 
they often suffer from instability and rheological problems. In 
particular, the particles tend to settle rapidly. Thus, although the 
slurries give better thermal conductivities, they are not practical. 
The nanofluid appears with development of nanoscience and nano- 
technology. The great potential of nanofluids in enhancing heat 
transfer means chances and challenges in thermal science and 
engineering. Since the concept of the nanofluid is newly proposed, 
so there are many questions which remain unclear and need to be 
solved [85,86]. 

Several important factors such as the type, shape and size of 
nano-particles, clustering and settling of nano-particles, and 
temperature of the nanofluids have not been studied widely. In 
addition, lack of agreement between experimental results and 
theoretical models is still the main issue of authors. 

Many authors reviewed the literature of nanofluids like 
[86-93] and reviewed and discussed the enhancement of con- 
vection heat transfer using nanofluids, experiments to synthesize 
the nanofluids, techniques of measuring the thermal conductivity 
and absolute viscosity of nanofluids, problems and challenges of 
nanofluids and the applications of nanofluids. In this review of 
literature, we focus on recent investigations which show new and 
precise empirical correlations for nanofluids properties with low 
deviations with related experimental results and new applica- 
tions for nanofluids particularly in micro-scale heat exchangers. 
4.1. Synthesis of nanofluids 

The synthesis of the nanofluids is the first key to enhance 
the thermal conductivity (k) of nanofluids. The convenient of 


nanofluid does not mean a simple liquid-solid mixture. Nano- 
fluids are produced by dispersing nanometer-sized solid particles 
into base-liquid such as water, ethylene glycol (EG), oil, etc. The 
major problem in the preparation of the nanofluid is the agglom- 
eration (clustering) of the nanoparticles. The nanotechnology 
opens wide opportunities to produce nanometer scale particle 
with different shapes. There are mainly two techniques for the 
preparation of nanofluids; the single (one)-step process and two- 
step process. 

The single-step process is direct evaporation approach. Pre- 
paration of nanofluid using this process is done by making and 
dispersing nanoparticles directly into the base fluid. This process 
is used to produce non-agglomerating copper nanoparticles 
which remain uniformly dispersed and stably suspended in EG. 
This process is favorable because it prevents the oxidation of the 
nanoparticle especially for nanofluids which have high-conduc- 
tivity metals such as copper. 

When the nanoparticles are made and then dispersed in the 
base fluid later it is known as two-step process. This process is 
extremely utilized to synthesize nanofluids considering the avail- 
able commercial nanopowders. The advantage of this process is 
that the inert gas condensation technique has been developed to 
economically produce nanoparticles. Whereas the disadvantage 
of this process is essentially when the nanoparticles tend to 
agglomerate before the nanoparticles dispersed in base fluid. So, 
one of the challenges to solve the agglomeration problem of the 
nanoparticles is the development of an innovative method to 
avoid the agglomeration which generally occurs in all nanopow- 
ders [94-96]. 

Liu et al. [97] prepared the nanofluids phase change materials 
(PCMs) wherein TiO, nanoparticles using the two-step procedure. 
The nanoparticles were uniformly dispersed in saturated BaCl, 
aqueous solution with pH 8. For stabilization, a supersonic 
oscillator was used to improve the dispersion of nanoparticles 
in fluid. They did not use any thickening agents to help suspend- 
ing and prevent phase separation. 

Zhu et al. [98] synthesized and developed stable CuO nano- 
fluids using a wet chemical method in a large-scale. They 
found that the concentration of copper acetate and reaction 
time affected the size and shape of clusters of primary nanopar- 
ticles. Nanofluids with different microstructures could be 
obtained through changing the synthesis parameters. The thermal 
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conductivity of CuO nanofluids increased with the increase of 
particle loading. It is expected that this method can be extended 
to synthesize other nanofluids. 

Sudhan and Meenakshi [99] developed the one-pot method 
for preparing silver nanofluids by reducing silver nitrate using 
sodium hypophosphite as reducing agent and ethylene glycol as 
base fluid by means of conventional heating. This is advantageous 
over other conventional methods; economical and one-step. The 
non-agglomerated and stably suspended silver nanofluids are 
obtained in a short time and the synthesized silver nanofluid 
has superior thermal conductivity. Hence, the silver nanofluid 
resulted, as they mentioned, may be used as an effective coolant 
in automobile industry in place of the conventional fluids that are 
currently in use. 


4.2. Mechanisms of nanofluids 


Keblinski et al. [100] explored four possible explanations for 
anomalous increase in the thermal conductivity with the decrease 
in the grain size. These mechanisms are Brownian motion of the 
particles, molecular level layering of the liquid at the liquid/ 
particle interface, the nature of heat transport in the nanoparti- 
cles, and the effects of nanoparticle clustering. They showed that 
the key factors in understanding thermal properties of nanofluids 
are the ballistic, rather than diffusive, nature of heat transport in 
the nanoparticles, combined with direct or fluid-mediated clus- 
tering effects that provide paths for rapid heat transport. 

Based on Green-Kubo linear response theory, Nie et al. [101] 
utilized the exact expression for the heat flux vector of the base 
fluid plus nanoparticle system to estimate the contribution of 
nanoparticle Brownian motion to thermal conductivity. They 
found that the contribution is too small to account for abnormally 
high reported values. They also found the possibility of convec- 
tion caused by Brownian particles is also unlike. They demon- 
strated a layer structure can form around the nanoparticles and 
the structure does not further induce fluid-fluid phase transition 
in the bulk fluid. For the models of an asymmetric hard sphere 
mixture representing the single spherical nanoparticles and a 
mixture of rods and hard spheres representing aggregates, both 
suspended in the fluid, they observed that for the very low 
volume fraction cases, the compressibility slightly changes. This 
shows that the speed of phonon transition is unchanged due to 
the addition of nanoparticles of any type. The results indicate 
that, besides the enhancement due to the high thermal conduc- 
tivity of nanoparticles themselves, fluid molecules make no 
evident contribution to the enhancement of thermal conductivity 
which is attributed to the presence of the nanoparticles at volume 
fractions less than 5%. 

Ma and Liu [102] demonstrated that having the highest 
conductivity, being electromagnetically drivable, the liquid metal 
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with low melting point is expected to be an idealistic base fluid 
for making super conductivity solution which may lead to the 
ultimate coolant in a wide variety of heat transfer enhancement 
area. 


4.3. Thermal conductivity measurement techniques 


The importance of thermal conductivity measurement is 
potential because k is responsible for heat transfer enhancement. 
Plenty of investigations have been done and developed in this 
area to get a reasonable accuracy. The temperature oscillation 
technique, the steady-state parallel-plate technique, and the 
transient hot wire method have been used to measure the 
thermal conductivity of nanofluids [103,104]. The last one has 
been utilized widely due to its high accuracy. The schematic 
diagram of the transient hot-wire set-up for measuring the 
thermal conductivity of nanofluids is shown in Fig. 18. 

Because of the electrical conductivity for all nanofluids, the hot 
wire should be surrounded by such material as epoxy adhesive 
which has an excellent electrical insulation to avoid applying the 
hot wire directly to the nanofluid [106]. There is no significant 
negative effect of the insulation layer on the performance of this 
method. The principle of this method is using a thin metallic wire 
embedded in the liquid to measure its thermal conductivity. This 
wire is employed as a heat source and a temperature sensor. The 
wire is heated by applying a power supply through it. The 
principle idea of this method is measuring the temperature/time 
response of the wire to an abrupt electrical pulse. A higher 
conductivity of surrounding nanofluid, a lower temperature will 
be raised in the wire. The advantage of this method is faster than 
the other techniques, and it has a simple design compared to 
other techniques. 

Xie et al. [103] measured thermal conductivity of nanosized 
SiC (26 nm) suspensions using the same method. Their experi- 
mental results illustrated that the thermal conductivities of the 
suspensions containing a small amount of solid particles are 
higher than those of the base liquids. For low volume fraction, 
they found that thermal conductivities increase almost-linearly 
with the volume fraction of the particles. The k increased by about 
15.8% at 4.2vol% of particles. For different systems containing the 
same kind of particles, the thermal conductivity ratios are 
independent of the base liquids, but the absolute value of the 
thermal conductivity of the suspension is proportional to that of 
the base liquids. Wang and Xu [107] measured the thermal 
conductivity of nanofluids consisted of Al2O0; (28 nm) and CuO 
(23 nm) suspended with distilled water as a base fluid by a 
steady-state parallel method. They demonstrated that the pre- 
dicted thermal conductivities of nanoparticle-fluid mixtures are 
much lower than their measured data. They indicated that there 
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Fig. 18. (a) Schematic diagram of the transient hot-wire apparatus for measuring the thermal conductivity of nanofluids [105] and (b) schematic diagram of transient hot 


wire instrument [54]. 
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is a deficiency in models when they employed for nanoparticle- 
fluid mixtures. 

Asseal et al. [104] obtained a 34% of enhancement of the 
thermal conductivity for 0.6vol% of carbon multi-walled nano- 
tubes (C-MWNTs) suspended in water using the transient hot- 
wire method. Hu et al. [108] measured thermal conductivity of 
AlN-ethanol nanofluids by a hot-disk method from 0.5% to 4% 
volume fraction at temperatures of 273.15 K and 297.15 K. The 
hot-disk thermal constant analyzer system is shown in Fig. 19. 

They reported that 20% increase in thermal conductivity of 
ethanol with the addition of 4vol% at 273.15 K and a strong 
temperature dependence of the thermal conductivity. Li et al. 
[109] measured the thermal conductivity of Al203 (47 nm)/water 
nanofluid with 0.5-6 vol% using the transient hot wire and the 
steady-state cut-bar method at room temperature. They con- 
cluded that the enhancement in thermal conductivity of the 
nanofluid exists and is independent of the measurement techni- 
que. Chen et al. [110] applied small angle X-ray scattering (SAXS) 
to the characterization of SiO, nanoparticles (10-30 nm) uni- 
formly dispersed in a water-base fluid. They observed that 
thermal conductivity for 16 vol% of nanoparticles shows a linear 
increase with increasing the average particle size. They suggested 
that thermal resistance at the solid-liquid interface may be 
critical to the thermal conductivity of nanofluids with small 
particle sizes. 

A transient Hot Wire Technique for measuring the thermal 
conductivity in liquids is implemented in an in-house made 
computer controlled experimental set-up was presented by 
Codreanu et al. [111]. The aim of this technique is that the 
thermodynamic state of the material under test is not changed 
during the measuring process due to the small heat disturbance 
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Fig. 19. Hot-disk thermal constant analyzer system, (A) computer, (B) computa- 
tion device, (C) Keithley 2000, (D) Keithley 2400, (E) Hot-disk bridge, (F) sensor, 
and (G) samples [108]. 
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and thermal properties can be assigned to the equilibrium 
temperature attained prior to the disturbance. They reported that 
the Quick Start Development System from Analog Devices offers a 
very convenient, rapid, low cost, highly sensitive and accurate 
solution for implementation of the Hot Wire Technique. The high 
resolution of the converter ensures a good precision of measure- 
ments. Fig. 20 illustrates the experimental set-up and the testing 
cell configuration. 

Table 1 shows the experimental results of the enhanced 
thermal conductivity of nanofluids measured by different techni- 
que methods. 


4.4. Thermal conductivity enhancement 


In 1990, Bonnecaze and Brandy [121] developed a general 
method to predict the effective thermal conductivity of an 
infinite, statistically homogeneous suspension of particles in an 
arbitrary (ordered or disordered) configuration. The method 
follows closely that of ‘stokesian dynamics’ and captures both 
far-field and near-field particle interactions accurately with no 
convergence difficulties. The results are found to be in excellent 
agreement with the accepted values. The effective conductivities 
of spherical particles in cubic arrays are calculated for particle to 
matrix conductivity ratios of infinity 10 and 0.01. Xuan and Li [94] 
observed that the volume fraction, shape, dimensions and proper- 
ties of the nanoparticles affect the thermal conductivity of 
nanofluids. They revealed that the nanofluid shows great 
enhancement in the heat transfer process. The ky,/k; of the Cu- 
H20 nanofluids remarkably increases with the volume fraction of 
ultra-fine particles. This increase is found to be varied from 1.24 
to 1.78 times when the volume fraction of particles increases 
from 2.5% to 7.5%. 

Suspensions of micron-sized silicon-carbide (SiC) particles 
with various aspect ratio ranged from 1 to 10 distributions are 
prepared and measured by Cherkasova and Shan [122]. They 
inferred that the thermal conductivity results are above the 
Maxwell prediction when anisotropic particle are used. They 
noticed, when the volume fraction remains constant, the thermal 
conductivity enhancement goes from 16.5% up to 39.5% as the 
volume-averaged aspect ratio increases from 4.8 to 9.6. They 
observed that any possible contamination of the suspension by 
the grinding material itself should only increase the conductivity 
with grinding time which is opposite to the observed trend. 
Kondaraju et al. [123] investigated the effective thermal conduc- 
tivity of Cu (100 nm)/water nanofluid and Al,03 (80 nm)/water 
nanofluid by developing the Eulerian-Lagrangian based direct 
numerical simulations (DNS) model. Their results exhibit an 
increase in the thermal conductivity of nanofluids with the 
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Fig. 20. (a) Basic configuration of the experimental set-up measuring the thermal conductivity of liquids and (b) the testing cell configuration [111]. 
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Table 1 


Illustration of the techniques used for measuring k. 
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Nanofluid/Ref. Diameter/length 
TiO2/BaClz aqueous [97] 20 nm 
SiC/water [103] 26 nm 

600 nm 


C-MWNTs/water 
C-DWNTs/water [104] 


Dour=40 nm, Din =15 nm 
Dour=5 nm, Djp»=2.5 nM 


1=10 um 
Al203/water 28 nm 
CuO/water [107] 23 nm 
ALN/ethanol [108] 20nm 
Al203/water [109] 47nm 
Si02/water 12nm 
Al203/water [112] 30nm 
Al,03/water 45 nm 
Al,03/EG [113] 
ALN/ethanol [114] 20nm 
Cu20/water 200.5 nm 
CuSO,/water [115] 
CuO/water 29 nm 
Al203/water [116] 36 nm, 47 nm 
Al,03/water 47 nm 
CuO/water [117] 
Silver-water [118] 100-500 nm 
Hybrid sphere of alumina/iron oxide CNTs [119] 070 nm, 24m 
AlzoCugo/EG [120] 20-40 nm 


Volume fraction 


Method used to measure k 


Max. enhancement in k 


0.167-1.13% THW 15.65% 
4.2% Transient hot wire 20% 
0.6% THW 34% 
38% 
2.4-8% Steady-state 40% 
Parallel-plate 
0.5-4% Hot-disk 20% 
0.5, 2, 4, 6% THW kruw < Ksteaay-state at 4 vol% 
Steady-state cut-bar krw = Ksteady-state at 2.6 vol% 
0.45-4% 30 
0.5-1.5% 
4% 30 13.3% 
0.5%-4% Hot-disk 20% 
0.02-0.05 THW 24% 
mol/L 
3.3-9.32% THW 16% (20-40 °C) 
3.14-9% 16% (20-40 °C) 
0.8% Transient hot wire 6.52% 
24.6%, and 
2-4 times at 30-60 °C 
1000 ppm Transient hot wire 10% 
5000 ppm 16% 
10 000 ppm 18% 
0.2% 3w 21% 
0.2-2 vol% Modified thermal comparator Twofold 


increase of volume fraction. They proposed that the temperature 
two-way coupling term should be considered for predicting the 
thermal conductivity of nanofluids. The results also show larger 
thermal conductivity for nanofluids where the coagulation of 
particles is not considered. The coagulation of nanoparticles in 
fluids is inevitable and thus numerical models should consider it 
for an accurate prediction of thermal conductivity. 

Yu and Choi [124] developed the Hamilton-Crosser model for 
suspensions of nonspherical particles to include the effect of a 
solid/liquid interface. The solid/liquid interface is described as a 
confocal ellipsoid with a solid particle. Their new model for the 
three-phase suspensions is mathematically expressed of the 
equivalent thermal conductivity and equivalent volume fraction 
of anisontropic complex ellipsoids as well as empirical shape 
factor which is expressed as follows respectively: 


i kp—ks 
kyj = {1 | TOD REA ow 


y (a2 +t\(b? +t +t) 
a abc f= 


J rf 2) 


2e(t)[1—e2(t)]'/° 


a e(t),/1—e2(t)+arcsin e(t) 


This model correctly predicts the magnitude of the thermal 
conductivity of nanotube-in-oil nanofluids. But at present, this 
model is not able to predict the nonlinear behavior of the 
nanofluid thermal conductivity. Shukla and Dhir [125] presented 
numerical results obtained from Molecular Dynamics Simulations 
of a solid-liquid system comprising of Lennard-Jones atoms used 
to study the liquid layering on solid nanoparticles. They found 
that close to the solid surface the liquid atoms from ordered 
layers which display higher thermal conductivity compared to the 
bulk liquid. They also presented a model for thermal conductivity 
of nanofluids based on the theory of Brownian motion of a free 


(3) 


particle as follows: 


F 11 bkg’T 
A particle = liquid + 3272 ath 


where u is the viscosity of the liquid. This relationship agrees 
qualitatively with the experimental results in that contribution of 
nanoparticles to the increase in thermal conductivity rises with 
an increase in particle volume fraction and fluid temperature. 

Kumar et al. [126] derived a model of thermal conductivity 
based on the statistical mechanics in the first term and Brownian 
motion based on critical size and optimum effective volume 
fraction in the second term as follows: 


keff kp +(N—1)kp—(n—1) e(Key w Cae(T—To) 
kp |  kp+(n-1)kf+&e(kf—kp) Unrkeré 


where %=0(1+h/r,)°, h is the liquid layer thickness and re is the 
critical particle size. The lesser the particle size, the higher the 
liquid layer thickness. 

Putnam et al. [127] claimed that it is not observed an 
anomalous enhancement of the thermal conductivity that has 
been reported in previous studies of nanofluids. The largest 
increase in thermal conductivity that they noticed is 1.3+0.8% 
for nanoparticles having 4nm in diameter and volume fraction 
<1 for Au nanoparticles suspended in ethanol. Murshed et al. 
[128] conducted a combined experimental and theoretical study 
on the effective thermal conductivity of titanium oxide and 
aluminum oxide nanoparticles suspended in DIW and EG. They 
emphasized that the effective thermal conductivity is substan- 
tially higher than the values of the base fluids. They proposed 
models for nanofluids which show a good agreement with the 
experimental results and give better predictions for the effective 
thermal conductivity. They found that the effective thermal 
conductivity of nanofluids increases with the nanoparticle 
concentration. 

Wang et al. [129] proposed a fractal model for predicting the 
effective thermal conductivity of liquid with dilute suspension of 
nonmetallic nanoparticles. This model predicts well the trend for 


(4) 


(5) 
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variation of the effective thermal conductivity with dilute sus- 
pension of nanoparticles, and fits successfully with their experi- 
mental data for CuO (50 nm)/DIW nanofluid with volume fraction 
less than 0.5%. Their calculated results show that the predictive 
calculation of effective thermal conductivity is complicated. 
Emami et al. [130] derived a general simple equation for thermal 
conductivity of all suspensions, including nanofluids, based on 
distance between particles. The Brownian motion approach 
should be used to calculate this distance in the nanofluids. The 
interactions between nanoparticles that were neglected in the 
derivation of Brownian motion relations may yield an incorrect 
calculated distance. An adjustable parameter was used to correct 
this distance as follows: 


a _ 71/2 
Abraw wri 2ct—-344e 2%] —dp (6) 


Han et al. [119] synthesized a new type of complex nanopar- 
ticle-hybrid sphere (@70nm, 2 m)-CNT particle, consisting of 
numerous CNTs attached to an alumina/iron oxide sphere as 
shown in Fig. 21. These hybrid nanoparticles are dispersed to 
poly-alpha-olefin with sonication and a small amount of surfac- 
tant to form stable nanofluids. The thermal conductivity of the 
fluids was measured by a 3~m-wire method over a temperature 
range 10-90 °C. They states that for such hybrid nanoparticles, 
heat is expected to transport rapidly from one CNT to another 
through the center sphere and thus leading to less thermal 
contact resistance between CNTs when compared to simple CNTs 
dispersed in fluids. CNTs have an extremely high thermal con- 
ductivity, but thermal resistance between the CNTs and the fluid 
has limited their performance in nanofluids. Their experimental 
results refer that the effective thermal conductivity of the fluids is 
increased by about 21% at room temperature for particle concen- 
tration of 0.2 vol%. 

Corcione [131] presented empirical correlation for predicting 
the effective thermal conductivity of nanofluids, based on a high 
number of experimental data available in the literature. It is 
found that the kkr increases as the nanoparticle volume fraction 
and the temperature are increased, and the nanoparticle diameter 
is decreased as shown in Fig. 22. The wide range of the thermal 
conductivity equation of validity (i.e., 10-150 nm nanoparticle 
diameter, 0.002-0.09 volume fraction and 294-324 K tempera- 
ture), makes such equation useful by the engineering point of 
view. The model formulas of thermal conductivity of [131-134] 
are as follows, respectively: 


keffa _ kp + 2kp—2ep(kp—kp) 
kf Kp +2kf + p(kp—Kp) 


(7) 
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Fig. 21. TEM image of a sample hybrid sphere/CNT particle produced by the 
aerosol method. The central spheres are polydisperse with a broad size range from 
10 to 300 nm, having a geometric mean diameter of 70 nm [119]. 
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Fig. 22. Distributions of kem)/kr vs. p for Al203-H20, with d, and T as parameters 
[131]. 
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where lpr is the mean free path of water, B, is Boltzman constant. 
kem) in Eq. (7) is the effective thermal conductivity of the 
nanofluid according to Maxwell. 

Meibodi et al. [135] proposed a model of thermal conductivity 
of the suspension considering Brownian motion and interfacial 
layer as well as a new mechanism considering nanoparticles 
as liquid-like particles. Unlike Beck et al. [136], they pointed out 
that this model can be used for estimating of upper and lower 
limits of nanofluid thermal conductivity without any adjustable 
parameter. They concluded that the most important factor for 
thermal conductivity enhancement of nanofluids is distance 
between nanoparticles. At very low particle concentration, ther- 
mal conductivity can be calculated using Brownian approach. 
For micro-particles and/or high particle volume fractions can be 
calculated using effective diameter. They reported that Brownian 
motion and vibration of nanoparticle mechanisms can be con- 
sidered for thermal conductivity of nanofluids, but these mechan- 
isms display minor roles for enhancement of thermal conductivity 
in nanofluids. 

The self-consisted scheme is generalized to predict the ther- 
mal conductivity of nanofluids containing spherical nanoparticles 
with a conductive interface was done by Hadjov [137]. He 
assumed that a flux jump in the particle-fluid interface in the 
opposite to the assumption for temperature jump in the case of 
thermal barrier resistance. An upper and lower bound to the 
homogenized suspension thermal conductivity according to the 
particle packing is. The model takes into account the enhanced 
conductivity and thickness of this thin layer. The suspension 
thermal conductivity is dependent on the kind of particle packing. 
This packing depends on the particle surface state and size 
distribution. In the same area of research, Xue [138] presented a 
novel model of the effective thermal conductivity for nanofluids 
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considering the interface effect between the solid particles and 
the base fluid in nanofluids. The results of effective thermal 
conductivity of nanotube/oil nanofluid and Al,03/water nanofluid 
are in good agreement with the experimental data. Their model 
can interpret the anomalous enhancement of the effective ther- 
mal conductivity of nanotube/oil nanofluid and its nonlinearity 
with nanotube loading. His formula of thermal conductivity is 


ke—Kex 


1D) ke—km Lv i ke—kcy 
S akak I eee 


=0 
2ke +(1 —Bz x)(Key—ke) 
(12) 


where A=abc/[(a+t)(b+t)(c+t)], kj is the effective dielectric 
constant component, and Bp; is the depolarization factor compo- 
nent of the elliptical particle along the j-symmetrical axis. 

Koo and Kleinstreuer [139] theoretically found that the impact 
of nanoparticle Brownian motion is much more significant than 
the thermo-phoretic and osmo-phoretic motion effects. They 
observed that Brownian motion effect on thermal conductivity 
decreases with particle size. At very dilute nano-suspensions 
< 0.5%, particle interaction is negligible, while it may be profound 
for > 1%. The thermophoretic and osmo-phoretic motion effects 
were found to be independent of particle size, which is not the 
case for Brownian-motion included elevated thermal conductiv- 
ities of dilute nanofluids. Their formula of thermal conductivity 
due to osmo-phoresis is 


o1 atk  3ke 
~ 3m Ue ka+2ke 


kos PaCaVT (13) 


With viscous effects is 


o1 ok  3ke 
~ 3 ue ka+2ke 


kos (1 x 10°p,c,) VT (14) 
where k, and kg are the thermal conductivities of the continuous 
and discrete phases, respectively, VT is the temperature gradient 
in the spherical particle. 

As a continuous work, Raykar and Singh [140] proposed a 
theoretical model employing a differential effective medium with 
Brownian motion to predict effective thermal conductivity (ETC) 
of CNT nanofluids. They inferred that ETC is influenced signifi- 
cantly by Brownian motion and enhancement is higher in dilute 
nanofluids. They predict that in future it is expected that the 
terms for the third phase (surfactants), dispersibility factors, and 
dynamic interfacial resistance are also incorporated using proper 
mathematical descriptions to predict the ETC of nanofluids. The 
combined static and dynamic mechanisms based model for 
predicting the ETC of nanofluids is presented also by Murshed 
et al. [141]. The model shows a good agreement with the 
experimental data of several types of nanofluids and gives better 
predictions compared to the classical and recently models. This is 
because they took into account the effects of particle size, 
nanolayer, Brownian motion, and particle surface chemistry and 
interaction potential which are the static and dynamic mechan- 
isms responsible for the enhanced effective thermal conductivity 
of nanofluids. Their combined model of ETC of nanofluids is 


i fp O(kp—coks) [2y} —y3 +1] + (kp + 200k, )y} [$p77(@-1)+ 1] 
e =)" 73 (Kp + 2eak,)—(kp + Op) py [V3 +731] 


16 2ko +3k 20 * 


1 3KgT(1—-1.5y3 Pp Gr 
+ [pacat aa a ate ey 


where d, is the distance between two particles, cp is the specific 
heat, Kg is the Boltzmann’s constant, r is the particle radius, and 4 
is the viscosity. For smaller sized nanoparticles and low volume 


34? 9A? Kept2kp  3.A4 
fargas T mE . 


fractions, dynamic mechanisms such as particle Brownian motion, 
particle interactions and surface chemistry are significant in 
enhancing the thermal conductivity of nanofluids. However, the 
major contributions to the enhanced thermal conductivity of 
nanofluids arise from static mechanisms. It can be inferred that 
the thermal conductivity of nanofluids is enhanced due to both 
static and dynamic mechanisms. 

Papari et al. [142] presented a neural network method to 
estimate thermal conductivity of nanofluids consisting of 
MWCNTs suspended in oil («-olfin), decene (DE), DW, EG, and 
single-walled carbon nanotubes (SWCNTs) in epoxy and poly 
methylmethacrylate (PMMA). Their predicted thermal conductiv- 
ities values are in good agreement with the other theoretical 
models as well as experimental results. They mentioned that 
examining the neural network approach to predict other trans- 
port properties such as viscosity of nanofluids remains for future 
work. 

Wang et al. [143] derived a new structural model of ETC for 
heterogeneous materials with multiple continuous phases by 
using mathematical deduction, a thermal field method, and an 
average field approximation. This new model has a distinctive 
structure which is substantially different from five conventional 
fundamental structural models (Series, Parallel, two forms of 
Maxwell-Eucken, Effective Medium Theory). The new model 
provides more structural models using the combinatory method 
proposed by literature. The model also provides narrower bounds 
of the ETC within the Hashin-Shtrikman bounds for heteroge- 
neous materials where the physical structure can be character- 
ized. The thermal conductivity of the CNTs composite consists the 
combined effect of CNT length, diameter, concentration, interface 
and matrix as 


91- ke—Km Lf ke—K$; kek, O 
2ke+km ` |ke+0.14(k§3—Kke)  2ke+ (1/2) (KS, ke) 
(16) 
where 
Cc ke Cc ke 
kn = FORJ 3 THOR ke/L) a7 


The results of Koroteeva et al. [144] demonstrated that the 
numerical technique based on the extended Maxwell methodol- 
ogy works well for materials with periodic and random arrange- 
ment of spherical cavities. The analysis indicates that, for the 
face-centered cubic (FCC) and body-centered cubic (BCC) array of 
spherical pores, the Maxwell approximation provides relatively 
accurate estimate for the ETC up to large porosities. For the 
simple cubic (SC) array, the Maxwell methodology accurately 
predicts the keg only for the porosities that are less than 0.3. The 
SC, BCC, and FCC schematic diagram are shown in Fig. 23. 

Li and Peterson [145] analyzed and modeled the mixing effect 
of the base fluid in the immediate vicinity of the nanoparticles 
caused by the Brownian motion. The simulation results indicate 


Fig. 23. (a) SC, (b) BCC, and (c) FCC [144]. 
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Fig. 24. (a) The velocity field of single nanoparticles simulation, (b) the comparison for simulated temperature field vertical plane (left, the temperature of X-Z plane; right, 


the temperature of X-Y plane) [146]. 


that this mixing effect can have a significant influence on the ETC 
of nanofluids as shown in Fig. 24. 

Koo and Kleinstreuer [147] modeled a new equation to calculate 
the thermal conductivity for the nanofluid including the Brownian 
motion. Their formula depends upon the particle volume fraction, 
particle size, particle material and temperature. Their formula is 


keff = Kstatic F KBrownian a 8) 
where 

(knp + are 
Kstatie =K 19 
K static f (knp +2kf)+ (kp—Knp) ( ) 


kT 
KBrownian =5x 10*BobpsCps Ipa TO (20) 
np 


While f in Eq. (20) is a correlation that depends on the type of 
nanoparticles and shown in Table 2. 


f(T, p) = —6.04 +0.4705)T +(1722.3ġ—134.63) 
for 1% <  <4% and 300 K < T < 325 K (21) 


Vajjha and Das [148] developed the model of Koo and 
Kleinstreuer [147] by deriving new model for f and f(T, ¢) from 
their set of experimental data from three nanofluids. Their new 
model of f is shown in Table 3. 

And f(T, p) is 


f(T,p) = (0.028217 +3.917 x 10°)(7-) +(—3.0669 


x 10776—3.91123 x 1077) (22) 
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Yurong et al. [149] correlated an equation for the thermal 
conductivity of TiO2/water nanofluid as follows: 


knf = ki(125.6267 +4.82+1.0) (23) 


4.5. Thermal stability 


Kuznetsov [150] proposed a novel type of a nanofluid that 
contains both nanoparticles and motile (oxytactic) microorgan- 
isms to obtain enhanced mass transfer, micro-scale mixing, and 
anticipated improved stability of the nanofluid. 

Sharma et al. [118] synthesized a silver nitrate-EG and poly 
nanofluids with different concentrations of silver nanofluid ran- 
ged from 1000 to 10,000 ppm. The results indicate a 10,000 ppm 
silver nanofluid exhibits rapid increase in the particle size with 
the passage of time as shown in Fig. 25. After 30 days of 
preparation, the thermal conductivity of 1000 and 5000 ppm 
silver nanofluids is observed to decrease slightly from 10% and 
15% to 9% and 14%, respectively. In addition, the thermal con- 
ductivity of 10,000 ppm nanofluid was decreased from 18% to 14% 
after 30 days. It is noted that the silver particles were aggregated 
in early stage of preparation (up to 15 days), which leads to the 
increase in the size of silver particles. However, no significant 
change was observed after 15 days, which indicates the stability 
of silver nanofluids. The thickness of dispersion stabilizer sur- 
rounded each particle is reduced which might be the reason of 
increasing aggregation of silver particles as the time lapsed. 

Sundar and Sharma [117] presented experimental results of 
heat retaining time of Al203 and CuO nanoparticles suspended in 


Table 2 
Correlations of f for different nanoparticles [147]. 


Type of particle B Remarks 

Au-citrate, Ag-citrate and CuO 0.0137(100@)- 9-229 p<1% 

cuo 0.0011(100¢)~°7272 b>1% 

Al,03 0.0017(100ġ) 70284 p>1% 
Table 3 


Correlations of f for different nanoparticles [148]. 


Type of particle f Concentration (%) Temperature (K) 


Al,O3 8.4407(100ġ) 19734 1<@<10 298 <T<363 
ZnO 8.4407(100)- 197304 1<@>7 298 <T<363 
Cuo 9.881(100¢)-°°45 1<ġ>6 298 <T < 363 


PSA average particle size (nm) © 


1 day 


7days 15days 30 days 
Time (day) 
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water with 0.8% nanoparticle volume fraction. Their measure- 
ments include measuring 50 ml in flask and heat up. The water, 
AlO; nanofluid, and CuO nanofluid have taken 17, 25 and 
27.4 min, respectively to reach 35°C. This analysis implies by 
adding some nano-size metallic oxide particles to the base fluid, 
these particles absorb heat and release heat to the adjacent liquid 
layers. But for a homogeneous fluid this will not be possible, 
because the entire liquid layer cool equally. Karthikeyan et al. 
[151] examined the behavior of thermal conductivity of CuO 
(8nm) suspended in water and EG base fluid with 1 vol% 
nanoparticle concentration. They highlighted that the time 
dependent thermal conductivity in water based CuO nanofluid 
shows that the thermal conductivity decreases with elapsed time 
due to clustering of nanoparticles with time, as confirmed 
microscopically. 


4.6. Viscosity 


Viscosity as a physical property of fluids is significantly 
influenced by the volume fraction of nanoparticles and tempera- 
ture of nanofluids. Viscosity is dramatically affected by Re number 
and friction factor of flowing fluid. Luciu et al. [152] found that 
the absolute viscosity of alumina/water nanofluid with volume 
fraction of 0.0-4.0% is reversely proportional to the temperature 
of the nanofluids as shown in Fig. 26. Whereas this viscosity is 
increasing as the nanoparticles concentration increases. 

Einstein [153] established the following equation to calculate 
the viscosity: 


Leff = Hy +2.5¢) (24) 
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Fig. 26. Variation of viscosity [152]. 
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Fig. 25. (a) Average particle size of silver nanoparticles as a function of lapsed time and (b) thermal conductivity of silver nanoparticles for different concentrations of 


silver [118]. 
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The correlations of viscosity of [131,132,134] are as follows, 
respectively: 


Leff 1 


= 25 
u 1-3487(dp/d) 7p ® (25) 
Leff = Hof F Happ (26) 
where {Happ is the apparent viscosity defined as 
p,Vpde 


where C depends on the base fluid viscosity and mean diameter of 
the nanoparticles, 6 depends on the mean diameter and volume 
fractions of the nanoparticles and vg is the Brownian velocity of 
the nanoparticles that depends on temperature, diameter and 
density of particles [154]. 


Herp = Ustatic + MBrownian (28) 
where 
Ustatic = My /0 -03 (29) 
kT 
a) 104 =———— fT, 30 
HBrownian x Pops 2DapRog $) (30) 


Sommers and Yerkes [155] displayed that the viscosity of 
aluminum oxide (10 nm) dispersed in propanol exhibits non- 
linear behavior with concentration with the viscosity increasing 
sharply for concentrations greater than 1.0%. Torii [156] found 
experimentally that the viscosity of nano-diamond fluid increases 
with increasing the nanoparticles concentration. He observed that 
the experimental data is higher than the theoretical values over a 
wide range of volume fraction. While he studied three types of 
nanofluids which are alumina, copper oxide and diamond nano- 
particles dispersed in conventional fluid of water [157]. He 
pointed out that the measured viscosities of nanofluids are much 
higher than that of predicted values using Batchelor equation 
which is 


Mit, 69 eae 4 (31) 
Uy 

The viscosity of nanofluids increases with an increase in the 
volume fraction and this trend is different for three nanoparticles. 
He stated that the effective volume fraction including cluster in a 
nanofluid becomes higher than that in the ideal suspension fluid 
in which each particle is independently and homogenously 
dispersed in a fluid. He has seen that the viscosity of diamond 
nanofluid is the highest because diamond particles are strongly 
aggregated and cluster restricts a large amount of pure water. On 
the contrary, the viscosity of CuO nanofluid is relatively low 
because the CuO nanoparticles had already aggregated in the 
state of the powder. 

Tavman et al. [112] illustrated that the effective viscosity of 
SiOz/water and Al,03/water has a similar behavior as water with 
the increase in temperature which ranged from 20°C to 50°C. 
Garg et al. [158] showed that the viscosity increase is almost 
fourfold of that predicted by the Einstein law of viscosity. 

Due to the large increase in viscosity than the increase in 
thermal conductivity, nanofluids would be poorer coolants as 
compared to base fluids if they are used in the existing heat 
transfer hardware. Schmidt et al. [159] found that the measure- 
ments of longitudinal viscosity and the corresponding values of 
the particle size of alumina suspended in polyalphaolefin (PAO) 
are consistent with a picture of non-clustered particles in weakly 
shear-thinning viscous oligomeric oil. 

Duangthongsuk and Wongwises [160] carried out experiments 
to investigate the effect of the thermophysical properties model 


on the predicted values of the convective heat transfer coefficient 
of TiO2-water nanofluid flowing in a horizontal double-tube 
counter flow heat exchanger under turbulent flow conditions. 
They showed that the Brinkman equation Eq. (31) and Einstein 
equation Eq. (32) give similar results of viscosity, whereas the 
Wang et al. equation Eq. (33) gives a higher predicted viscosity 
than the preceding formulas. 


1 
Unf = a-p (32) 
Ung =(14 Py 63) 
Mn =(1+7.3ġ +1230) Hw (34) 


Chopkar et al. [120] presented experimentally that CaCO3- 
distilled water nanofluid with different nanoparticle concentra- 
tions is advantageous to reducing aggregation of primary nano- 
particles. The effective viscosities and effective thermal 
conductivities of the nanofluid are related to the aggregates of 
nanoparticles and can be well predicted by the modified Krieger 
and Dougherty formula and the modified Hamilton and Crosser 
model, respectively. Whereas the modified Krieger and Dougherty 
formula for viscosity is 


Nesp /Ny = (1-@4/P_q) 7 
Pa = P(a/ay? (35) 


where a, and a are the ratio of aggregates and primary particles, 
respectively, D is the fractal index, has a typical value of 1.8 for 
nanofluids. 

Murshed et al. [128] conducted the combined experimental and 
theoretical study on the effective thermal conductivity and viscosity 
of titanium oxide and aluminum oxide nanoparticles within 1-5 vol% 
dispersed in DIW and EG. The viscosity of nanofluids is observed to be 
substantially higher than the values of the base fluids. Both thermal 
conductivity and viscosity of nanofluids increase with the nanopar- 
ticle concentration. Corcione [132] presented empirical correlation for 
predicting the dynamic viscosity of nanofluids, based on a high 
number of previous experimental data. 


Hepp _ 1 
My = 1-34.87(dp/dp) 97 p18 


(36) 


where d; is given by 


6M 1/3 
s- (it) en 


In which M is the molecular weight of the base fluid, N is the 
Avogadro number, and pfo is the mass density of the base fluid 
calculated at temperature T,=293 K. Corcione [161] and Cianfrini 
et al. [162] developed Eq. (36) by multiplying it by a factor of 0.1. 

In addition, the (Hu) increases as the nanoparticle volume 
fraction is increased, and the nanoparticle diameter is decreased, 
being practically independent of temperature for a wide regions of 
the dynamic viscosity equation of validity (ie, 25-200nm in 
diameter, 0.0001-0.071 volume fraction and 293-323 K tempera- 
ture), makes such equation useful by the engineering point of view. 
Sundar et al. [163] experimentally and theoretically examined the 
effect of temperature variation on several flow characters such as 
density, specific heat and kinematic viscosity. They demonstrated 
that the density and viscosity of nanofluids decreases with the 
temperature for various concentrations of nanoparticles. While the 
specific heat increases with the temperature. 

Yurong et al. [149] presented their experimental data of 
viscosity of the TiOz/water nanofluid at about 22°C by the 
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following correlation: 


Hn = Mpp(199.217 +4.62$ + 1.0 (38) 


While Nguyen et al. [164] proposed a formula for 4 vol% of CuO 
nanoparticles suspended in water base fluid as follows: 


Unf = Hwater(2.1275—0.0215T +0.0002T) (39) 


And for Al,03/water with dp=47nm and volume fraction 
ranges from 0% to 13% 


Mnf _ 9.904480) (40) 


Log 

where T is the temperature in °C. Vajjha [165] developed a general 
correlation for viscosity of Al203, CuO and SiO2 nanoparticles 
suspended in (60%EG/40%W) as follows: 


Eng =A,e%?) (41) 


Log 

where A; and A2 are constants. This equation works with 
specifications of nanoparticles diameters and volume fraction 
ranges. 


4.7. Temperature dependence of thermal conductivity enhancement 


Most measurements of the former investigators have been 
done at room temperature without attention to that the tem- 
perature has an important effect on the particle movement which 
tremendously affects on the thermal conductivity. Das et al. [166] 
were the first one to experimentally examine the thermal con- 
ductivity dependence on the temperature oscillation. They carried 
out their tests on water as a base fluid with particles of Al,03 or 
CuO as suspension material. Their results bring out that a 
dramatic enhancement in thermal conductivity takes place with 
temperature. They obtained an enhancement of thermal conduc- 
tivity up to 2% with 1% particles of Al,O3 at room temperature 
21 °C, while a 10.8% enhancement was obtained at 51 °C. While a 
24.3% enhancement was seen with a 4% volume fraction for the 
same previous range of temperature. They attributed this aug- 
mentation occurred in thermal conductivity of nanofluids to the 
increase occurred in the Brownian motion of particles. They 
stated that nanofluids can be made more attractive as cooling 
fluid for devices with high energy density where the cooling fluid 
is likely to work at a temperature higher than the room tempera- 
ture. The temperature dependence of thermal conductivity 
enhancement of water-Al,03 nanofluids can be seen in Fig. 27. 
This trend of thermal conductivity was also seen by Beck et al. 
[167] for alumina dispersed in either water or EG. 

Minsta et al. [116] carried out experimental tests to study the 
relationship between the thermal conductivity and temperature 
oscillation. Particles of CuO with in 29 nm as well as 36 nm and 
47 nm of Al203 particles is dispersed in water to produce nano- 
fluids. The temperature range was from 20 °C to 40 °C. They found 
that the ETC of nanofluids increases with temperature for the 
three types of nanofluids. It is observed that the linear tendency 
fits of the measurements has slope more than 10°. They obtained 
16% enhancement in thermal conductivity for each type of 
nanofluids when temperature ranged between 20°C and 40°C. 
The effect of temperature oscillation on thermal conductivity 
depends on the value of particles concentration. In comparison 
with distilled water, they found that an addition of nanoparticles 
gives a better enhancement with temperature for low volume 
fraction of particles. 

Sundar and Sharma [168] stated that the ambient temperature 
thermal conductivity property of nanofluid is not sufficient for 
estimating the heat transfer coefficient. Thermal conductivity at 
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Fig. 27. Temperature dependence of thermal conductivity enhancement for 
water-Al,03 nanofluids [167]. 
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Fig. 28. Thermal conductivity of Al,O3 nanofluid at different temperature [168]. 


different temperature is very important. They found that thermal 
conductivity of alumina/water and copper oxide/water nanofluids 
are relatively proportional with temperature shown in Fig. 28. 
Depend upon the increase of temperature, thermal conductivity 
of nanofluid increases with volume fraction. 

The same results were obtained by Izadi et al. [169] for the 
alumina—water nanofluid with temperature range from 300 K to 
325K. Liu et al. [97] obtained a maximum enhancement in 
thermal conductivity for TiO2 nanoparticles dispersed in satu- 
rated BaCl aqueous solution with pH 8 which was 15.65% with 
1.13 vol% in higher temperature of 15 °C. They stated that the 
concentration of nano-particles and temperature are two impor- 
tant factors in improving thermal transport properties of 
nanofluids. 

Li et al. [170] studied the effects of several parameters on the 
thermal conductivity of CuO (29 nm)/DW and Al203 (36 nm)/DW 
with 2-10 vol% of nanoparticles and nanofluid temperature ran- 
ged from 27.5°C to 34.7°C. Their results indicate that the 
temperature has a significant impact on the ETC. At temperature 
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of 34°C and 6 vol%, they observed an enhancement in thermal 
conductivity by about 1.52 times for CuO/DW nanofluid greater 
than that of DW. At the same previous condition of temperature 
and with 10 vol%, they remarked 1.3 times of thermal conductiv- 
ity enhancement for Al,03/DW nanofluid higher than that of the 
base liquid. Wei et al. [115] showed that the thermal conductivity 
of nanofluid exhibits a strong sensitivity to the temperature. Beck 
et al. [167] measured the thermal conductivity of alumina 
nanoparticles dispersed in water, EG, and (water+EG) mixture 
nanofluids over a temperature range from 296 K to 420 K. They 
predicted the thermal conductivity and the effects of particle size, 
volume fraction, temperature and changes in base fluid using a 
simple model. They showed that the thermal conductivity- 
temperature behavior follows closely that of the base fluid. The 
maximum in thermal conductivity occurs near 380 K, which is 
close to the maximum in thermal conductivity-temperature 
behavior for pure EG or (water+EG) mixture. These results 
confirm the hypothesis that the temperature dependence of the 
effective thermal conductivity of nanofluids is due mostly to that 
of the base fluid. 

Bianco et al. [171] studied numerically the turbulent forced 
convection flow of Al,0O3/water nanofluid in a circular tube 
subjected to a constant and uniform temperature at the wall. To 
simulate the nanofluid convection, two-phase mixture model is 
considered taking in account appropriate thermophysical proper- 
ties. They assumed that the particles are spherical with a 38 nm in 
diameter. They obtained a good agreement with the experimental 
correlation proposed in literature. They stated that the accuracy 
of the model could be improved by using a better description of 
nanofluid thermophysical properties taking into account the 
temperature dependence. Murshed et al. [128] conducted a 
combined experimental and theoretical study on the effective 
thermal conductivity of titanium oxide and aluminum oxide 
nanoparticles suspended in DIW and EG. They found that the 
thermal conductivity of nanofluids is strongly dependent on 
temperature. The same results are observed by Yu et al. [172] 
for the copper (5-10 nm)/EG nanofluids with 0.5 vol% and at 
50°C. They have seen that the thermal conductivity strongly 
depends on the temperature of fluid, and the enhancement ratios 
of the thermal conductivity of nanofluid increase along with the 
increasing temperatures. Brownian motions of Cu nanoparticles 
would play the key role in determining the effect of the tempera- 
ture on thermal conductivity enhancement of nanofluids. 

Temperature dependence of thermal conductivity of Alz03/EG 
nanofluid was studied by Beck et al. [136] for a wide range of 
temperature from 298 K to 411K. The k of the nanofluid was 
measured by liquid metal transient hot wire apparatus. They 
confirmed that additional temperature contributions inherent in 
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Brownian motion models are not necessary to describe the 
temperature dependence of the thermal conductivity of nano- 
fluids. They showed that the effect of mass or volume fraction of 
nanoparticles on the thermal conductivity of nanofluids can be 
correlated using the Hamilton and Crosser or Yu and Choi models 
with one adjustable parameter (the shape factor in the Hamilton 
and Crosser or the ordered liquid layer thickness in the Yu and 
Choi model). 

Sundar and Sharma [117] presented the increase of thermal 
conductivity with temperatures for nanofluids with water as base 
fluid and particles of Al203 or CuO as suspension material. It is 
observed that a twofold to fourfold increase in thermal conduc- 
tivity enhancement of nanofluids can take place over a tempera- 
ture range of 30-60 °C. It was observed that nanofluids containing 
smaller CuO particles show more enhancement of conductivity 
with temperature. Sankar et al. [173] proposed a theoretical 
approach based on molecular dynamic modeling to estimate the 
thermal conductivity enhancement of platinum-water nanofluid. 
The results distinctly show that the enhancement of thermal 
conductivity increases as the volume fraction of nanoparticles 
increase from 1% to 7%. The enhancement is found to be more for 
higher temperatures. The temperature effect of thermal conduc- 
tivity enhancement in Al203/water and CuO/water nanofluids was 
also presented by Sundar and Sharma [174]. They observed a 
dramatic increase in the enhancement of thermal conductivity 
occurs with temperature as shown in Fig. 29. 


4.8. Effect of volume fraction 


The previous studies demonstrated that there is a linear 
relationship between the thermal conductivity of nanofluids and 
the concentration of nanoparticles. But the increase in concentra- 
tion of nanoparticles leads to undesirable consequents such as 
clogging, sedimentation, and pressure drop. Kumar et al. [175] 
found numerically a linear dependence of thermal conductivity 
enhancement on particle concentration of Au nanoparticle (e < 1) 
with 4nm in diameter. Liu et al. [97] found that the thermal 
conductivity of Tiz0/saturated BaClz aqueous nanofluids increases 
approximately linearly with the volume fraction of nanoparticles. 
A similar trend of thermal conductivity they observed with the 
changing the nanofluid temperature in different volume fraction 
of nanoparticles. This tendency of thermal conductivity with the 
volume fraction of nanoparticles is also observed by Jain et al. 
[176] which are similar to that of previous literature. Wei et al. 
[115] reported that the variation in the CuSo, molar concentra- 
tion leads to a change in nanofluid thermal conductivity through 
changing the nanofluid microstructure. A linear relationship 
between the ratio of thermal conductivity of alumina/water 
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nanofluid and the concentration of nanoparticles was also proved 
by Izadi et al. [169]. Hu et al. [114] revealed that 20% enhance- 
ment in thermal conductivity of alumina/ethanol nanofluid can be 
obtained when 4.0% of nanoparticle is dispersed in the base liquid 
at 273.15 K. The same enhancement ratio in thermal conductivity 
was obtained by Hu et al. [108] where 4.0% volume fraction of AIN 
nanoparticles with 20 nm in diameter is dispersed in ethanol at 
273.15 K. 

Xie et al. [103] emphasized that the thermal conductivities of the 
SiC suspensions containing a small amount of solid particles are 
significantly higher than those of the base liquids. At low volume 
fraction range, the thermal conductivity enhancing ratios increase 
almost-linearly with the volume fraction of solid particles. The 
experimental results of Sunder and Sharma [168] showed that water 
containing a small amount of alumina or copper nanoparticles does 
not exceed 0.8% provides a substantially higher thermal conductivity 
than the same liquids without nanoparticles. Tavman et al. [112] have 
found experimentally that the ETC of alumina and silica nanoparticles 
dispersed in water increases as the volume fraction increases but not 
anomalously as indicated in literature. They claimed that this increase 
is independent on temperature. Labonté et al. [177] found numeri- 
cally that for variable properties model, the enhancement in heat 
transfer is 20% where 1% particles concentration of alumina dispersed 
in a base liquid of water and 40% with 4% particles concentration. 
While for constant properties model, they obtained 5% and 20% for 
the same previous particles concentrations, respectively. 

An increase of 1.52 times in the effective thermal conductivity 
of CuO nanoparticles with 6 vol% dispersed in DW was obtained 
by Li and Peterson [170]. This enhancement of the ETC of alumina 
nanoparticles was screened while 10% volume fraction is dis- 
persed with water. The enhancements were at temperature of 
34°C. 

Many parameters respect to the volume fraction were exam- 
ined by Sundar et al. [163]. These parameters are the density, 
specific heat, and kinematic viscosity of alumina/water nanofluids 
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with a volume fraction up to 0.8% of alumina particles. They 
screened that when metallic nanoparticles is added to the base 
fluid such as water those particles act as suspensions and the 
density of the water increases. This is because the density of the 
nanofluids is directly proportional to the mass of the fluid. They 
showed that the specific heat is reversely proportional to the 
volume fraction of nanoparticles. Their results revealed that the 
viscosity of the nanofluid increases with the increase of the 
nanofluid concentration. These results are shown in Fig. 30. 

Das et al. [166] showed that the thermal conductivity of 
alumina-water and CuO-water nanofluid is directly proportional 
with the concentration of the nanoparticles at different tempera- 
tures. But this is stronger with alumna than with copper oxide. 
The same tendency of thermal conductivity with nanoparticle 
concentration was also seen by Jang and Choi [178] at volume 
fraction less than one percent for CNTs and less than 0.05% for an 
oxide and metallic nanoparticles. Sommers and Yerkes [155] 
showed that the enhancement of thermal conductivity of an 
alumina nanoparticles with 1 vol% dispersed in propanol was 
found to be slightly higher than that of base liquid for 
1800 < Re < 2800 and at room temperature. They attributed that 
to an earlier transition to turbulent flow. For 3% volume fraction, 
they observed that the thermal performance is deteriorated with 
respect to the baseline case. 

The effective thermal conductivity of Ar-Cu nanofluid in shear 
field was calculated by equilibrium molecular dynamics (EMD) 
simulation using Green-Kubo formula was achieved by Sun et al. 
[179]. The relative thermal conductivity increases linearly as the 
flow shear rate increases for the enhanced “micro-convection” 
effect yielded by the rotation of the nanoparticle. The increase in 
the ETC at lower volume fraction with respect to the shear rate is 
larger than at higher volume fraction which is related to the 
slowed down rotation speed for the higher volume fraction of 
nanoparticle. The results exhibit that the shear flow did not affect 
the fluid structure and orderly fluid structure is not the reason for 
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Fig. 30. (a) The density, (b) the specific heat, and (c) the Kinematic viscosity of nanofluid at different concentrations [163]. 
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the enhancement in the thermal conductivity of shearing nano- 
fluid. The effective thermal conductivity obtained from the con- 
ventional correlation: 


= 1 +BcPep 
to 
where 2e is the effective thermal conductivity of the flowing suspen- 
sion, Ag is the thermal conductivity of the suspension at static state, c 
is the volume fraction of the particles, B and m are constants, and Pep 
is the Peclet number of the particle expressed by 


(42) 


_ ed? 
Op 
where e is the velocity gradient, d is the particle diameter, and « is 
the thermal diffusivity of the suspended fluid. It is obviously seen that 
the volume fraction effect is more explicit for the flowing suspensions 
containing micro-sized particles is significantly lower than their 
numerical results. Therefore, they stated that the conventional 
correlation is not suitable when the sizes of the suspended particles 
are reduced to nanometer. 

Sundar and Sharma [117] reported that the thermal conduc- 
tivity of Al,O3 and CuO nanofluids at 0.8 vol% increases by 6.52% 
and 24.6% when compared to base fluid of water. The enhance- 
ment is considerably increased for nanofluids with Al2O3 as well. 
Sharma et al. [118] synthesized nanofluids consisting silver 
nitrate (precursor), EG (reducing agent), and poly (acrylamide- 
co-acrylicacid) (dispersion stabilizer). The results illustrate an 
increase in thermal conductivity of silver nanofluids up to 10%, 
16%, and 18% as the amount of silver particles in nanofluid is 
1000, 5000, and 10,000 ppm, respectively. The effect of the 
nanoparticle concentration of Al,03/water and CuO/water nano- 
fluids was also implemented by Sundar and Sharma [174]. They 
highlighted that the thermal conductivity is directly proportional 
with the volume fraction as shown in Fig. 31. They presented that 
thermal conductivity of previous nanofluids for 0.8 vol% increase 
by 6.52% and 16.00%, respectively compared to the base fluid. 
They showed that nanofluids containing smaller amount of CuO 
particles show more enhancement of thermal conductivity com- 
pared to AlO; particles. 

Karthikeyan et al. [151] reported that the increase in the 
enhancement of thermal conductivity of CuO (8 nm) suspended 
in water and EG base fluid with 1 vol% was observed to be 31.6% 
and 54%, respectively. They stated that the thermal conductivity 
nonlinearly increases with the volume fraction of nanoparticles. 


Pep (43) 


4.9. Effect of particle size 


The decrease in the nanoparticles size means an increase in the 
surface area of the nanoparticles exposed to the fluid. There is a 
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significant effect of the nanoparticle size on the thermal con- 
ductivity of nanofluids as shown in the literature. Kumar et al. 
[175] displayed the increases in the enhancement of the effective 
conductivity with the reciprocal of the nanoparticles radius as 
shown in the following equation: 


_ Keff 
km(1—€)rp km 


q kpêrm 
~ —kmAm (dT/dx) ~ 


q“ (44) 
where the subscript m refers to the liquid and p refers to the 
nanoparticles. The results of Li and Peterson [170] reveal that the 
thermal conductivity of CuO (29 nm)/water is much higher than 
that of Al203 (36 nm)/water because the copper oxide nanopar- 
ticles is smaller than that of Al2O3 nanoparticles. It is greatly 
confirmed by the new model of Jang and Choi [178]. Their new 
model is 


Kopp = Ker —f) + Bkparticief + C1 ns ker Rez, no Prf (45) 
where C,=18x 10°, f is the volume fraction, and kgr is the 
thermal conductivity of the base fluid molecules. Fig. 32 shows 
the nanoparticle diameter decreases as the effective thermal 
conductivity of nanofluids increases. 

As the particle size decreases, the Brownian motion of nano- 
particles is greater and then nonconvection becomes dominant. 
This consequent is identical to the results of Vasu et al. [180]. 
They analytically observed that the effective thermal conductivity 
of nanofluids increases as the volume fraction increases, the nano- 
fluid temperature increases and the nanoparticle size decreases as 
shown in Fig. 33. The same trend of nanofluid thermal conductivity 
was seen by Jain et al. [176]. The latter shows that the size effect is 
the most important effect in nanofluids, which traditional models 
fail to predict. 

The results obtained by Izadi et al. [169] show that the effective 
thermal conductivity of alumina/water nanofluids is also dramati- 
cally affected by the size of nanoparticles until 5 nm in diameter of 
nanoparticles. At larger than this size, they claimed that the effect 
of the nanoparticle size becomes nil. Beck et al. [181] performed 
experiments on alumina/water and alumina/ethylene glycol nano- 
fluids and they claimed that the thermal conductivity enhance- 
ment decreases as the particle size decreases below about 50 nm. 
The model made by Kumar et al. [126] showed that the Brownian 
motion contribution is higher at lesser particle size when particle 
volume fraction is held constant. As shown in Fig. 34, the changes 
in Brownian motion is in the range of 10-15 nm and this range 
is the critical size due to the considerable change in Brownian 
motion when the volume fraction and temperature are fixed. 
The Brownian motion is suppressed when the particle size and 
volume fraction are increased. This may be due to the fact that the 
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Fig. 31. (a) Thermal conductivity of AlO; nanofluid at ambient temperature and (b) thermal conductivity ratio of CuO nanofluid at ambient temperature [174]. 
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Fig. 32. Effect of diameter of nanoparticle on thermal conductivity of nanofluid [178]. 
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Fig. 33. The effect of (a) volume fraction, (b) temperature and, (c) diameter of nanoparticle on thermal conductivity of nanofluid [180]. 


Brownian motion velocity is inverse dependence with the particle 
diameter [126]. 

While Kondaraju et al. [182] studied the effect of multi-sized 
nanoparticle distribution on thermal conductivity of nanofluids 
ranged from 100 to 300 nm were considered while an initial 
diameter of 182.4 nm was considered for the single-sized nano- 
fluid case. Their simulated results of Al203/DIW nanofluid 
revealed an increase in the effective thermal conductivity for 
the nanofluids with multi-sized nanoparticles when compared to 
that of the nanofluids with the corresponding single-sized nano- 
particles. Sundar and Sharma [117] presented the increase of 
thermal conductivity with temperatures for nanofluids with 
water as base fluid and particles of Al2O3 or CuO as suspension 
material. However, the enhancement is considerably increased for 


nanofluids with Al,03 as well. The measurements indicate that 
particle size is an important parameter for the observed behavior. 
Sundar and Sharma [174] confirmed that the thermal conductiv- 
ities of AlpO3;water and CuO/water nanofluids increase as the 
particle size decreases. 

Lu and Fan [183] implied that the simplified dynamics simula- 
tion method is an effective method to forecast some thermal 
properties of nanofluids for volume fraction less than 8%. Their 
results give the effects of the volume fraction and the size of 
nanoparticles on the thermal conductivity and the viscosity of 
nanofluids. It is observed that the decrease in the nanoparticle 
size or increase in the volume fraction causes an increase in 
thermal conductivity and viscosity. For suitable volume fraction 
and size, increasing viscosity with improving heat transfer 
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capability is acceptable. For example, the largest enhancement of 
thermal conductivity is found to be 41% for AlO, (10 nm)/EG 
nanofluid with nanoparticle concentration of 5%. The effect of 
volume fraction on thermal conductivity is shown in Fig. 35. 

Karthikeyan et al. [151] reported that the increase in the 
enhancement of thermal conductivity of CuO (8 nm) suspended 
in water and EG base fluid with 1 vol% nanoparticle concentration 
has observed to be 31.6% and 54%, respectively. They attributed 
this dramatic increase in thermal conductivity to the finer particle 
size and monodispersity of nanoparticles. Saterlie et al. [184] 
obtained an enhancement of thermal conductivity up to 22% over 
water for the 0.55 vol% Cu nanofluids. When the particle loading 
was increased to 1.0 vol%, the nanofluids of oleic acid-prepared Cu 
powders settled and clogged the test setup. A maximum increase 
was 48% for the 1.0 vol% copper nanofluid from the cetyl 
trimethylammonium bromide (CTAB) prepared Cu powders. This 
enhancement in thermal conductivity is attributable to the 
excellent dispersion of the nanoparticles in the fluid. 


4.10. Effect of particle shape 


The thermal conductivity enhancement of nanofluids is com- 
pared with the geometrical shape of the nanoparticles as shown 
in Figs. 36 and 37. As shown in these figures, a large difference is 
existed when the geometrical shape is varied. The whiskers 
particles of SiC exhibits a large enhancement of thermal con- 
ductivity compared to the spherical particles as shown in Fig. 36. 
In addition, the CNTs particles provides good enhancement of 
thermal conductivity with respect to spherical geometry particles, 
while the hybrid particles of sphere/CNTs particles provide better 
augmentation in thermal conductivity as shown in Fig. 37. 
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Fig. 34. Variation of particle size with thermal conductivity ratio due to Brownian 
motion [126]. 
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4.11. Effect of particle material 


The particle material has significant influence on the thermal 
conductivity enhancement based on the conductivity of the 
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Fig. 36. Relative thermal conductivity enhancement in the suspensions containing 
SiC spherical particles (solid squares) and whiskers (solid circles) compared to the 
data calculated with the aid of Maxwell (dashed line) and EMT (solid line). Data by 
Xie et al. for spheres and cylinders are shown with hollow squares and circles 
correspondingly. EMT prediction for a=3.7 is shown with dotted line for 
comparison with the data on cylinders of Xie et al. [122]. 
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Fig. 35. Thermal conductivity as a function diameter of (a) Al,03/water and (b) Al,03/EG with 5 vol% [183]. 
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material of nanoparticle. The effect of the particle material 
has been studied by [94,123,137,156,185-187]. As shown in 
Fig. 38, the thermal conductivity of Al,03/EG, TiO2/EG, and 
Al203/H20, Cu/H20 and CuO/H20 is varied depending on the 
particle material. 


4.12. Effect of base fluid material 


The base fluid material plays an important role in the enhance- 
ment of the thermal conductivity of nanofluids [107,143,187]. 
Fig. 39a, shows that the alumina—water provides a better enhance- 
ment in thermal conductivity than the EG, engine oil, and pump 
fluid, respectively. The same tendency can be seen with the CuO 
nanoparticles suspended in EG and water as shown in Fig. 39b. 

The effect of the particle material and base fluid material as 
mentioned in the last two paragraphs can be concluded by the 
investigation done by Xie et al. [188]. They examined the 
enhancement of thermal conductivities of the DW, EG, glycerol, 
silicone oil, and the binary mixture of DW and EG as base 
fluids, and Al2O0; with different sizes, SiC with different shapes, 
MgO, ZnO, SiOz, Fe304, TiO}, diamond, and CNTs nanoparticles 
with different pretreatments were used as additives. The results 
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demonstrate and confirm that the thermal conductivity enhance- 
ments of nanofluids could be influenced by multi-faceted factors 
including the volume fraction of the dispersed nanoparticles, the 
tested temperature, the thermal conductivity of the base fluid, the 
size of the dispersed nanoparticles, the pretreatment process, and 
the additives of the fluids. In addition, the ky/ky decreases while 
the ky increases. 


4.13. Convection heat transfer 


An experimental study to investigate the increase in the 
thermal conductivity of the Al,03/DIW nanofluid under develop- 
ing convective boundary layer condition in tubes of 5mm in 
diameter was performed by Kolade et al. [146]. Their experiments 
were done in a hydrodynamically fully developed laminar tube 
flow in the range 500 < Re < 1600 with constant wall heat flux. In 
their study, the increase in effective thermal conductivity is 6% for 
2% volume fraction. For a suspension of MWNTs in silicone oil, the 
thermal conductivity is increased by 10% over that of the base 
fluid for 0.2% particle concentration. Based on the results of the 
carbon nanotube (CNTs) experiments, they suggested that the 
structure of the nanomaterial is an important factor in the 
augmentation of thermal conductivity of a liquid by the means 
of nanoparticles. 

Rostamani et al. [189] analyzed numerically the turbulent flow 
of mixture of CuO, Al,O3, and TiOz nanoparticles suspended in 
water with different nanoparticles concentrations flowing 
through a two-dimensional duct under constant heat flux condi- 
tion. All the thermophysical properties of nanofluids are tem- 
perature-dependent. The results imply that by increasing the 
volume fraction, the wall shear stress and heat transfer rates 
increase. For a constant volume fraction and Re number, the effect 
of CuO nanoparticles to enhance the Nu number is better than 
Al,03 and TiO nanoparticles.The effect of nanoparticle volume 
concentration on Nu number can be seen in Fig. 40. 


4.13.1. Forced convection heat transfer 

Lotfi et al. [190] studied numerically the forced convective of 
Al,03/water nanofluid flows in horizontal tubes. As a conse- 
quence, comparison of calculated results with experimental 
values shows that the mixture model is more precise than the 
single-phase and the Eulerian model. It is clearly shown that the 
rate of thermal enhancement decreases with the increase of 
nanoparticle concentration. Both single-phase and Eulerian 
approaches underestimate the Nu number in comparison with 
experimental results. Lai et al. [191] documented the results of 
heat transfer and pressure drop of alumina—water nanofluid 
flowing in a single 1.02 mm diameter stainless steel tube. The 
tube was electrically heated to provide a constant wall heat flux 
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Fig. 39. (a) Thermal conductivity as a function of volume fraction of Al2O0; powders in different fluids, (b) thermal conductivity as a function of volume fraction of CuO 


powders in ethylene glycol and water [107]. 
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Fig. 40. The influence of the nanoparticle volume concentration on the Nusselt number over a range of the Reynolds numbers [189]. 
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Fig. 41. The influence of Al,O3 nanoparticle volume concentration on the Nusselt 
number [192]. 


and the flow was under laminar flow condition. The heat transfer 
coefficient was measured in both developing and fully developed 
regions. Their experimental results illustrate that the heat transfer 
coefficient increases with the volume flow rate and nanoparticle 
volume fraction. In the developing region, the h enhancement 
decreases with increasing axial distance from the test section 
entrance. It is pronounced that the higher the volume fraction the 
longer the thermal entrance length. 

Sundar and Sharma [192] emphasized that Nu number is 
relatively proportional with Re number as shown in Fig. 41. This 
trend of Nu number with Re number is also confirmed by Xuan 
and Li [193]. 


Maiga et al. [194] investigated the laminar forced convection 
flow of yAl,03 suspended in the water and EG flows in a 
uniformly heated tube and a system of parallel coaxial and heated 
disks. It is pronounced that the inclusion of nanoparticles into the 
base fluids has produced an augmentation of the heat transfer 
coefficient that explicitly increases with an increase of the volume 
fraction. It is clearly shown that yAl,03-EG nanofluid provides a 
better heat transfer enhancement than yAl203-water. The 
enhancement in heat transfer increases with Re number. They 
have proposed correlations based on the numerical results for 
computing the averaged Nu number as a function of the govern- 
ing parameters Re and Pr: 


Nu = 0.086Re°°° Pr°> for constant wall heat flux (46) 


Nu = 0.28Re°?> Pr°36 for constant wall temperature (47) 


4.13.2. Natural convection heat transfer 

Polidori et al. [185] studied the natural convection heat 
transfer of alumina/water nanofluid up to 4% volume fraction. 
As a result, it was found that natural convection heat transfer is 
not solely characterized by the nanofluid effective thermal con- 
ductivity and that the sensitivity to the viscosity model used 
plays a key role in the heat transfer behavior. 

Jou and Tzeng [195] examined numerically the natural con- 
vection heat transfer enhancement utilizing nanofluids filled 
within a two-dimensional enclosure. The results indicate that 
increasing the buoyancy parameter and volume fraction of nano- 
fluids cause an increase in the average heat transfer coefficient. 
The numerical predictions illustrate the enhancement is critical 


5982 H.E. Ahmed et al. / Renewable and Sustainable Energy Reviews 16 (2012) 5951-5993 


0 
0 025 05 075 1 
X 


Fig. 42. Streamline and isotherm contours for different volume fraction (Gr=10°, Pr=6.2) [195]. (a) ¢=0, (b) ¢=0.1 and (c) ġ=0.2. 


for nanofluids than pure fluids. The stream line and isotherm 
contours for different volume fractions can be seen in Fig. 42. 

Putra et al. [196] studied the natural convection of nanofluids 
inside horizontal cylinder heated from one end and cooled from 
the other. They apparently observed in their experimental study 
the paradoxical behavior of heat transfer deterioration. The fluid 
shows characters distinction from that of common slurries and it 
appears to behave more homogeneously compared to common 
slurries. They documented that no evidence of long initial flow 
instability or stratification of concentration layers was found. 
Unlike conduction or forced convection, a systematic and definite 
deterioration in natural convective heat transfer was found to 
occur. The deterioration is dependent on particle density, con- 
centration as well as the aspect ratio of the cylinder. The reason 
for this effect is unclear. They stated that the role of particle-fluid 
slip and sedimentation seems to be important. Their study 
indicates the need for more investigations to understand the 
physical phenomenon completely before using these fluids in 
practical cooling applications. 


4.13.3. Boiling heat transfer 

Khandekar et al. [197] designed, fabricated and tested a two- 
phase gravity assisted thermosyphon with three nanofluids based 
on water with 1 vol% of Al203, CuO, and laponite clay. According 
to their results, it was observed that thermal performance 
deteriorates when nanofluids are used as working fluids. Max- 
imum deterioration was observed with laponite while minimum 
deterioration was for aluminum oxide particles based nanofluids. 
Increased thermal conductivity of the nanofluids does not appear 
to be affecting the nucleate pool boiling heat transfer coefficient 
although single-phase heat transfer improves. Boiling character- 
istics of nanofluids seems to be more effected by the physical 
interaction of nanoparticles with the nucleating cavities. 

Migration characteristics of nanoparticles in the pool boiling 
process of nanorefrigerant and nanorefrigerant-oil mixture are 
fundamental knowledge for the application of nanorefrigerants in 
refrigeration systems. Dinga et al. [198] experimentally studied 
and numerically simulated the migration characteristics of nano- 
particles in the nanorefrigerant-oil mixture as well as in the 
nanorefrigerant. Experimental results show that the migrated 
mass of nanoparticles in the pool boiling process of both nanor- 
efrigerant and nanorefrigerant-oil mixture increases with the 
increase of the original mass of nanoparticles and the mass of 


refrigerant. The migration ratio decreases with the increase of 
volume fraction of nanoparticles. The migrated mass of nanopar- 
ticles and migration ratio in the nanorefrigerant are larger than 
those in the nanorefrigerant-oil mixture. A numerical model, 
which qualitatively well predicted the migrated mass of nano- 
particles, was established, and the deviations between the model 
predictions and experimental data were in the range of 7.7-38.4%. 


4.14. Applications of nanofluids 


The demand to thermal systems has an ultra-high perfor- 
mance encouraged the researchers to focus their efforts to 
improve the energy efficient by using such fluids have high 
convective heat transfer coefficients as nanofluids. Nowadays, 
the nanofluids are being used in a wide range of industrial and 
commercial applications. A high performance of heat exchanger 
means a small size of the exchanger, light weight, low cost and 
diminishing in the fuel consumption in vehicle industries. 

Walsh et al. [199] demonstrated that -PIV measurements are 
feasible in nanofluids and can significantly aid in their under- 
standing, either by validating/negate many theories put forward 
to explain their enhanced transport properties or by providing the 
direction to develop more informed ones. 


4.14.1. Heat pipes 

Tsai et al. [200] studied the performance of nanofluids of 
aqueous gold nanoparticles solution used as a working fluid in 
heat pipe with different sizes and shapes of nanoparticles. Their 
results indicate that the thermal resistance of heat pipe is 37% for 
a nanofluid contains 0.2 ml of Na; citrate, 2.5 ml of Tannic acid 
and 3 ml of HAuCL, lower than that using DIW. They attributed 
this reduction in thermal resistance to that the suspended 
nanoparticles tend to bombard the vapor bubbles formed at the 
liquid-solid interface during their formation. So they reported 
it is expected that the nucleation size of vapor bubble is 
much smaller for fluid with suspended nanoparticles than that 
without them. 

Chen [201] examined the effectiveness of existence of silver 
nanoparticles with 35 nm in diameter suspended with pure water 
filled flat heat pipe (FHP). He found a small temperature differ- 
ence of FHP when a small amount of Ag nanoparticles added in 
the pure water. A 0.6 °C, 0.79 °C, and 1.24 °C temperature differ- 
ence is shown at 5 ppm of nanoparticles and 20-40 W was 
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utilized. These values are decreased when the nanoparticles 
concentration increases until 100 ppm as much as he tested. His 
results also illustrate that the thermal resistance of FHP filled 
with pure water is higher than nanofluid under 20-40 W. 

The efficiency of a two-phase closed thermosyphon using 
Al203 ( < 100 nm)/water nanofluid was improved by Noie et al. 
[202]. Different volume fractions of nanoparticles 1-3% are 
suspended within the two-phase closed thermosyphon (TPCT). 
The results inspire that nanofluids in all concentration studied 
show better thermal performance than pure water. The efficiency 
of the TPCT has improved up to 14.7%. It is distinctly seen that the 
temperature distribution on the TPCT are lower compared to pure 
water. Temperature differences between the evaporator and 
condenser sections with nanofluids are less than that pure water. 
On the other hand, thermal resistance of the TPCT when charged 
with nanofluids is less. This finding makes nanofluid attractive as 
working fluid in heat pipe and thermosyphon technology that 
further investigations are needed. 

Naphon et al. [203] studied the performance of heat pipe with 
de-ionic water, alcohol, and titanium with 21 nm in diameter 
suspended in alcohol. They observed a significant effect of the 
nanofluids on the enhancement of heat transfer due to higher 
heat capacity, higher thermal conductivity of working fluid and 
higher mixing fluctuation. They demonstrated that the heat pipe 
thermal efficiency increases with increasing nanoparticles con- 
centration. It was seen that 10.6% enhancement in heat pipe 
thermal efficiency is achieved while a 0.1 vol% of nanoparticles is 
used compared to that of the base fluid. When the volume 
fraction of nanoparticles exceeds 0.1%, the properties of the 
nanofluids seem to be a solid phase thereby lower evaporation 
rate of working fluid in the evaporator section. This in turn 
reduces the thermal efficiency of heat pipe. 
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Shafahi et al. [204,205] performed analytical models to study 
the thermal performance of a rectangular and disk-shaped heat 
pipes using Al,O3, CuO, and TiO2/water nanofluids with particle 
diameter ranged from 10 to 40 nm. They claimed that the liquid 
pressure drop decreases by using the nanoparticles within the 
water due to an increase in density. After reaching a critical 
concentration level, a reversal effect occurs resulting in an 
increase in the pressure drop as the increase in viscosity of the 
nanofluid. They screened a decrease in the velocity in both disk 
and rectangular shaped heat pipe when the concentration level 
increases. They showed that as the nanoparticle concentration 
increases the thermal conductivity of the nanofluid increases and 
then in turn the evaporator temperature decreases. They 
remarked that the smaller nanoparticles or larger concentration 
level the more pronounced the reduction in the temperature 
difference between the evaporator and the condenser. They have 
seen that the thermal resistance becomes smaller when a nano- 
fluid is utilized. It can be reduced up to 83% of its initial 
resistance. They stated that the CuO/water nanofluid provides 
30% and 20% reduction in the size of the disk and rectangular 
shaped heat pipes, respectively. The removal heat capacity from 
the disk-shaped heat pipe is shown in Fig. 43. 

The same previous nanofluids with a cylindrical heat pipe are 
tested by Shafahi et al. [205]. They found that the nanofluids 
enhance the heat pipe’s thermal performance by reducing the 
thermal resistance as the concentration increases or as the 
particle diameter decreases. They reported that smaller particles 
have pronounced effect on the temperature gradient along the 
heat pipe. Wei et al. [206] used an Ag/water nanofluid as a 
working fluid in a grooved circular heat pipe with average 
diameter of 10nm. Their test results show that the average 
decrease of 30-70% in the thermal resistance of the heat pipe 
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Fig. 43. Three dimensional representation of disk-shaped heat pipe maximum heat removal capability in terms of nanoparticles concentration levels and dimensionless 


thickness of the wick; (a) CuO, (b) Al203, and (c) TiOz [204]. 
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with nanofluid as compared with pure water. Tsai et al. [200] 
examined a gold/water nanofluid as a working fluid in a conven- 
tional circular heat pipe. They designed the heat pipe as a heat 
spreader for CPU in a notebook or a desktop PC. They showed that 
at the same charge volume, there is a significant reduction in 
thermal resistance of heat pipe with nanofluid as compared with 
DI water. They reported that the thermal resistance of a vertical 
meshed heat pipe varies with the size of gold nanoparticles. 


4.14.2. Electronics cooling 

With the fast improvement of computer performance and the 
CPU chips become so small, a dramatic heat has been generated in 
the chips need to escape using such effective working fluid as 
nanofluids. Nguyen et al. [207] investigated experimentally the 
heat transfer performance by using the Al203-water nanofluid for 
a liquid system that is utilized for cooling of electronic compo- 
nents. They have found that convective heat transfer coefficient 
increases as much as 23% with respect to that of the base fluid 
with 4.5% nanoparticle volume concentration. It was observed 
from their results that an augmentation of particle concentration 
has produced a clear decrease of the junction temperature 
between the heated component and the water-cooling block. 
The results exhibit that for particular mass flow rate of 0.04 kg/s, 
the enhancement of the thermal performance is about 10% and 23% 
for nanofluids with 2.2% and 4.5% particle concentrations, respec- 
tively when compared to the corresponding results of distilled 
water. 

Escher et al. [208] performed a systematic characterization of 
aqueous silica nanoparticle suspensions. They evaluated the 
potential of nanofluids for cooling of electronics by means of 
characterization of aqueous suspensions of silica nanoparticles for 
different volume fractions of 5-31%. They fabricated test-vehicles 
of microchannel heat sinks with three different channel widths. 
They found that an enlargement of the effective thermal mass 
which can be either achieved by increasing the heat capacity of 
the fluid or its density seems to be more promising to enhance the 
heat sink performance. 

Jung and Yoo [209] proposed a novel expression for the 
thermal conductivity of nanofluids which incorporates the kinetic 
theory to describe the contribution of the Brownian motion of the 
nanoparticles with more realistic definition of the mean free path. 
Their expression also considers the contribution of the interpar- 
ticle interaction due to the existence of the electrical double layer 
(EDL). They showed that this model is applied to Au/water 
nanofluids satisfactorily with respect to temperature, volume 
fraction and particle size. In the case of dense Al203/water, the 
effect of the interparticle interaction due to EDL on enhancing the 
thermal conductivity is more prominent than in the case of dilute 
Al203/water nanofluids. They outlined that this model shows that 
interparticle interaction due to EDL is the most responsible for the 
enhancement of thermal conductivity of nanofluids. 


4.14.3. Microchannels 

Tsai and Chein [210] performed an analytical study to examine 
the copper-water and carbon nanotube-water nanofluids as 
coolant in microchannel heat sink (MCHS). They found that the 
nanofluid reduces the temperature difference between the MCHS 
bottom wall and bulk nanofluid compared with that from pure 
fluid. This reduction of temperature is proportional to the particle 
volume fraction. This means a reduction will occur in conductive 
thermal resistance and an increase in convective thermal resis- 
tance when nanofluids are used due to the increase in viscosity 
and decrease in thermal capacity. They observed that the nano- 
fluid has a significant effect on the MCHS performance when the 
channel aspect ratio and porosity are low. Shokouhmand et al. 


[211] studied analytically the influence of Cu-water nanofluid 
with various volume fraction on the thermal performance of a 
silicon MCHS. They emphasized that nanofluids could enhance 
MCHS performance compared with that using pure water as a 
coolant. The enhancement is due to the increase in thermal 
conductivity of coolant and the nanoparticle thermal dispersion 
effect. They observed that there is no extra pressure drop 
recorded as long as the nanoparticle is small and particle volume 
fraction is low. 

Bhattacharya et al. [212] studied numerically the laminar 
forced convection heat transfer of Al203/H20 nanofluid flowing 
in a silicon MCHS. They have seen that using nanofluid improves 
MCHS performance by reducing fin conductive thermal resis- 
tance. The improvement of MCHS performance becomes more 
pronounced with increase of nanoparticles concentration. How- 
ever, they reported that there is an increase in pumping power by 
using nanofluids in MCHS which somewhat offsets the beneficial 
effect of using nanofluid. They revealed that the use the same 
nanoparticle concentration could give better synergy between 
velocity and temperature gradient vector at lower Reynolds 
number. 

Ho et al. [213] conducted experiments to investigate the forced 
convective cooling performance of a copper MCHS with alumina/ 
water nanofluid as the coolant in laminar regime. They remarked 
a significant high average heat transfer coefficient and markedly 
lower thermal resistance and wall temperature at high pumping 
power. Despite they marked increase in dynamic viscosity of 
nanofluid, the friction factor for the nanofluid was found slightly 
increased. It is observed an increase by about 70% in the average 
heat transfer coefficient with 1% nanoparticle concentration at 
Re=1676 compared with that with water. Lee and Mudawar 
[214] studied experimentally the effect of the same nanofluid 
implemented by Ho et al. [213] on the thermal performance of 
microchannel. Their results show a significant enhancement in 
thermal conductivity when nanofluid is used especially for 
laminar regime. High heat transfer coefficients were achieved 
mostly in the entrance region of microchannels. The enhance- 
ment was weaker in the fully developed region proving that 
nanoparticles have a remarked effect on thermal boundary layer 
development. Higher concentrations produced greater sensitivity 
to heat flux. They reported that despite this enhancement, the 
overall cooling effectiveness of nanoparticles was quite slight 
because of the large axial temperature rise associated with 
the decreased specific heat for the nanofluid compared to the 
base fluid. 

Lee et al. [215] studied the effective convection coefficient and 
viscosity in microtubes (D=0.5 mm) along with stationary ther- 
mal conductivity measurements for (Al,03, ZnO, and CuO) nano- 
fluids and carbon nanotubes (CNT). They obtained an effective 
convection coefficient increase of 5% for 3 vol% of Al,03/DIW 
nanofluid, 13.3% for 4 vol% of CuO/DIW nanofluid, and 11.6% for 
0.2 vol% for CNT/DIW nanofluid. 

Mohammed et al. [216] performed a 3-D simulation to study 
the effect of using Al203-H20 nanofluids on heat transfer and 
fluid flow characteristics in rectangular microchannel heat sink 
(MCHS). Their results revealed that when the volume fraction of 
nanoparticles is increased under the extreme heat flux, both the 
heat transfer coefficient and wall shear stress are increased while 
the thermal resistance of the MCHS is decreased. Nanofluid with 
volume fraction of 5% could not be able to enhance the heat 
transfer or performing almost the same result as Ovol% they 
found that a slight increase in the pressure drop in the MCHS 
compared with the pure water-cooled MCHS. While Gunnase- 
garan et al. [217] studied numerically in 3-D the effect of 
geometrical parameters which are rectangular, trapezoidal, 
and triangular on water flow and heat transfer characteristics in 
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microchannels. It is found that better uniformities in heat transfer 
coefficient and temperature can be obtained in the heat sinks 
having the smallest hydraulic diameter. It is also inferred that the 
heat sink having the smallest hydraulic diameter has better 
performance in terms of pressure drop and friction factor among 
other heat sinks studied. Mohammed et al. [218] performed a 
simulation to investigate the laminar flow and heat transfer 
characteristics of trapezoidal MCHS using various types of base 
nanofluids and various MCHS substrate materials on MCHS 
performance. they displayed that the best uniformities in heat 
transfer coefficient and temperature among the four mixture 
flows can be obtained using glycerin base nanofluid followed by 
oil-base nanofluid, EG-base nanofluid, and water-base nanofluid 
heat sinks. However, the heat transfer performance of water-base 
nanofluid can be greatly enhanced in steel made substrate 
heat sink. 


4.14.4. Heat exchangers 

Macro-compact heat exchangers have been widely used in 
various application thermal fluid systems including automotive 
thermal fluid systems, radiators for engine cooling systems, 
evaporators and condensers for HVAC systems, oil coolers and 
inter coolers. They are typical examples that can be found in 
vehicles have an important role in its weight and design of its 
front end module which has a strong impact on the car aero- 
dynamic behavior [187]. 

Sundar et al. [219] studied experimentally the effect of the 
alumina/water nanofluid flowing through a copper tube subjected 
to constant heat flux boundary condition in the turbulent regime 
with and without tape insert. From their results, it was seen that 
further enhancement in heat transfer with twisted tape are 
achieved when compared to that in plain tubes using nanofluid 
at the same flow rate. Farajollahi et al. [186] measured the heat 
transfer characteristics of y-Al2O3/water and TiO2/water nano- 
fluids in a shell and tube heat exchanger under turbulent flow 
condition. Based on the results, adding the nanoparticles to the 
base fluid causes a significant enhancement of heat transfer 
characteristics. For both nanofluids, two different optimum 
nanoparticle concentrations are found. Comparison of heat trans- 
fer behavior of two nanofluids indicates that at a certain 
Peclet number, heat transfer characteristics of TiO2/water nano- 
fluid at its optimum nanoparticle concentration are higher than 
those of y-Al2,03/water nanofluid. Whereas the latter nanofluid 
possesses a better heat transfer behavior at higher nanoparticle 
concentrations. 

Hwang et al. [220] examined experimentally the heat transfer 
and pressure drop characteristics of alumina/water nanofluids 
flowing through a uniformly heated circular tube in the fully 
developed laminar regime. From their results, it was shown that 
the results for nanofluid friction factor show a good agreement 
with analytical predictions from Darcy’s equation for single-phase 
flow. However, the convective heat transfer coefficient of the 
nanofluids increases up to 8% at a concentration of 0.3 vol% 
compared with that of pure water and this enhancement cannot 
be predicted by the Shah equation. Furthermore, it is observed 
that the convective heat transfer enhancement exceeds the 
thermal conductivity enhancement. 

Vasu et al. [187] performed a theoretical analysis with e-NTU 
rating method using also Al,03/water nanofluid as coolant on 
automobile flat tube plain fin compact heat exchanger. The results 
indicate that nanofluids possess higher heat transfer character- 
istics than conventional coolant water. The overall heat transfer 
coefficient is higher for nanofluids than water and increases with 
an increase in the nanoparticles concentration. They highlighted 
that with the increase of the volume fraction the cooling capacity 


in moderate manner and pressure drop increase with coolant 
inlet temperature. But cooling capacity is high when compared 
with pure water. It can be seen that in Fig. 44. 

Chopkar et al. [120] prepared a nanofluid containing a small 
amount ( < 1 vol%) of nanocrystalline metallic particle (AlzgCu39) 
as a stable colloidal dispersion in EG. Their study showed that the 
thermal conductivity is significantly greater 1.2-2 times than that 
of the base fluid as well as that predicted by the models. They 
pointed out that this enhancement of thermal conductivity is a 
function of size, volume fraction and thermal property of the solid 
suspension. A simple quenching experiment demonstrates that 
the efficiency of the nanofluid is better than that of the base fluid 
particularly at lower temperature. Hence, this nanofluid could be 
useful in microthermal heat transfer applications including in 
MEMS and automobiles. 

Mapa and Mazhar [221] performed an experimental study to 
investigate the effect of the nanofluid on heat exchanger thermal 
performance. The nanofluid implemented in their experiments 
was 99% pure copper oxide pre-dispersed in water with an 
average particle size of 29nm. They found that this nanofluid 
enhances the heat transfer rate. An increase in heat transfer rate is 
observed at any given flow rate. For example, at a mass flow rate 
of 0.005 kg/s, 5.5% increase in heat transfer rate is observed. As 
the concentration of nanoparticles in the fluid increases, a further 
increase in the heat transfer rate is observed. The presence of 
nanoparticles in base fluid reduces the thermal boundary layer 
thickness. Khoddamrezaee et al. [222] studied the influence of 
Al203-EG nanofluid and EG fluid which flow in a rectangular 
arrangement of tubes in a shell and tubes heat exchanger. They 
indicated that the use of alumina/EG nanofluid increases the heat 
transfer coefficient ratio by about 3.25 on the first column and 
3.16 on the second column. 

Advancement in the electronics industry led to the develop- 
ment of micro-scale heat transfer devices which offered high heat 
transfer coefficient in a compact size. The amount of research done 
in this particular field is fairly new and limited. Most studies done 
on microchannel devices and nanofluids recently have reported 
enhanced heat transfer capabilities and results that challenge 
traditional theories and limitations on heat transfer devices and 
fluids [223]. Mohammed et al. [224] studied numerically the 
effects of using various types of nanofluids on heat transfer and 
fluid flow characteristics in a square shaped of an aluminum 
microchannel heat exchanger (MCHE). They reported that the 
thermal properties and performance of the heat exchanger were 
enhanced using nanofluids while a slight increase in pressure drop 
was observed. They added that increasing the Re number causes an 
increase in the pumping power and a decrease in the effectiveness. 
While Mohammed et al. [225] found that an enhanced perfor- 
mance with the usage of nanofluids and slight penalty in pressure 
drop. The increase in Re number caused an increase in the heat 
transfer rate and a decrease in the overall bulk temperature of the 
cold fluid. The increase in nanoparticle volume fraction results 
better performance at the expense of increased pressure drop. 


4.14.5. Solar collector 

Heat transfer enhancement in solar devices is one of the key 
issues of energy saving and compact designs. The use of the 
additives is a technique applied to enhance the heat transfer 
performance of base fluids [226]. Natarajan and Sathish [226] 
concluded that the thermal conductivity enhancement depends 
on the volume fraction of the suspended particles and thermal 
conductivities of the particles and base fluids. Their experiments 
proved that the nanofluid is more effective than the conventional 
fluids. If these fluids are used as a heat transport medium, they 
increase the efficiency of the traditional solar water heater. 
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Fig. 44. Air and coolant flow influence on the thermal and fluid dynamic performance of compact heat exchanger [187]. 


Shahi et al. [227] numerically simulated the steady natural 
convictive flow and heat transfer in a single-ended tube with non- 
uniform heat input and with different inclination angle. This 
sample serves as a simplified model of the single-ended evac- 
uated solar tube of a water-in-glass evacuated tube solar water 
heater. The sealed end of the tube is assumed to be adiabatic and 
the opening end is assumed to be subjected to copper—water 
nanofluid with particle concentration up to 0.05%. This study 
indicates that the maximum overall mean Nu number is obtained 
at 35° inclination angle, while the maximum output mass flow 
rate is increasing with the inclination angle. It is also obtained 
that both the Nu number and maximum output mass flow rate 
increasing function of solid concentration, but the presence of 
nanoparticles is more effective at the smaller inclination angles. 
Shin and Banerjee [228] tested two types of nanoparticles powder 
fine grained powders and coarse powders. They demonstrated 
that the electron microscopy of the samples reveal that the 
variation in thermophysical properties of the different samples 
was caused by the amount of the agglomeration of the nanopar- 
ticles as well as the formation of compressed phase. They showed 
that these nanomaterials can enable significantly reduction in the 
cost of solar thermal power. 


4.15. Main problems of nanofluids 


Despite the effectiveness of nanoparticles in enhancing the 
single-phase heat transfer coefficient by increasing the nanofluid 
thermal conductivity, especially for laminar flow, but there are 
some disadvantages of nanofluids such as; possible damage to 


flow loop by erosion, and inability to sustain flow boiling [214]. 
Other disadvantages can be summarized as follows: 


4.15.1. Pressure drop 

In 2002, Qiang and Yimin [229] carried out an experimental 
investigation to study the convection heat transfer and flow 
structure of the Cu-water nanofluid flows through a tube in 
laminar and turbulent regime. When the nanoparticles suspended 
in the water at lowervolume fraction the friction factor has not 
change. Therefore, they stated that nanofluid with so small 
nanoparticle concentration behaves like pure water as shown in 
Fig. 45 [229]. Compared with water, no significant augmentation 
in pressure drop for the nanofluid is found in all runs of the 
experiments [193]. 

The pressure drop of the CNTs suspended in distilled water 
and flows in a horizontal tube was also studied by Ko et al. [230]. 
They found that under the laminar flow conditions the friction 
factor of CNTs nanofluids stabilized by adding surfactant is much 
larger than that of CNTs nanofluids prepared by acid treatment, 
and both nanofluids show larger friction factors than distilled 
water. Under the turbulent flow conditions, the friction factor of 
both nanofluids becomes similar to that of the base fluids as the 
flow rate increases as can be seen in Fig. 46. 

Sundar and Sharma [192] investigated numerically the turbu- 
lent fully developed flow heat transfer coefficient and friction 
factor of AlzO3/water and Al203/EG nanofluids in circular tube. 
They outlined that a 1.42 times of friction factor was achieved for 
a 4 vol% of nanoparticles and Re number of 30,000 over the base 
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Fig. 45. (a) The friction factors of nanofluids for laminar flow and (b) the friction factors of nanofluids for turbulent flow [229]. 
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Fig. 46. Friction factor as a function of dimensionless flow rate [230]. 
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Fig. 47. Comparison of pressure drop of Al,03 nanofluid in turbulent region [192]. 


fluid. They observed that pressure loss increases with increase in 
the volume fraction of the nanoparticles as shown in Fig. 47. 
Bontemps et al. [231] performed experiments to investigate 
the behavior of SiO2/water nanofluids inside a circular tube with 
imposed wall temperature. It was observed that high nanoparticle 
concentrations lead to surface modification probably due to 
particle deposition. This effect has to be confirmed for lower 
nanoparticle concentrations. In the turbulent regime, the pressure 
drop rapidly increased as a function of time to reach a maximum 


value. So they reported that Blasius law was not applicable 
anymore. They have seen an associated pressure drop to the 
presence of the nanoparticle in the base fluid. They observed the 
ageing effect just after loading the loop and after 40 h. Torii [157] 
revealed that the pressure drop of the nanofluids is slightly 
increased compared with that of the pure water, because an 
increase in the friction loss is caused by suspension of nanopar- 
ticles in the pure fluid. Lee and Mudawar [214] reported that the 
increasing in pumping power due to great pressure drop, and 
long-term fluid settling and potential clogging of flow passages 
are disadvantages of nanofluids. 


4.15.2. High viscosity 

The absolute viscosity of alumina/water nanofluid is increas- 
ing as the nanoparticles concentration increases and this trend is 
different for three nanoparticles [154-157]. Due to the large 
increase in viscosity than the increase in thermal conductivity, 
nanofluids would be poorer coolants as compared to base fluids if 
they are used in the existing heat transfer hardware [158]. The 
effective viscosities and effective thermal conductivities of the 
nanofluid are related to the aggregates of nanoparticles [232]. 


4.15.3. High cost of nanoparticles 

Because of the ultra-fine techniques (nanotechnology) and 
high accurate equipment used to synthesize and produce nan- 
ometer size, the nanofluid has a high cost. So, the high price of the 
nanofluid is one of the major barriers to apply nanofluids in a 
wide range of applications up to date. Lee and Mudawar [214] 
reported the cost of nanofluids is one of the disadvantages. 
These factors restrain the researchers to enlarge the range of 
applications. 
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4.15.4. Low specific heat 

As shown in literature, the specific heat of nanofluid is 
reversely proportional with the nanoparticles concentration 
[155,160,219]. If the concentration is more than the amount heat 
absorbed by the nanofluid [163]. Higher specific heat capacity of 
nanofluid, yields higher heat removal. An increased axial rise in 
wall temperature due to degraded specific heat is reported in 
negative side of nanofluids by Lee and Mudawar [160]. 


5. Conclusion 


The heat transfer and flow structure characteristics are studied 
with the presence of vortex generators (VG). Different types of 
VGs are displayed. The VGs are embedded, attached, punched or 
mounted in laminar or turbulent boundary layer were investi- 
gated. The impact of VG on convective heat transfer coefficient 
and friction factor when VGs are present in rectangular channel, 
triangular duct, flat-plate, fin-tube heat sink, and plate heat 
exchanger were reviewed. Several arrangements for single, row, 
and two-dimensional array of VGs through the flow direction are 
summarized. The main conclusions can be briefly drawn as 
follows: 


e Longitudinal vortices have as significant effective than the 
transverse vortices on heat transfer enhancement. 

e Pointing down rectangular winglet vortex generators yield 
higher heat transfer rate than pointing up rectangular winglet 
vortex generators. 

e Common flow-down VG gives better heat transfer character- 
istics than the common flow-up VG when a pair of DWVG is 
embedded in turbulent boundary layer of rectangular channel. 

e Winglets provide better performance in heat transfer than 
wings. Delta-winglet pair performs slightly better than a 
rectangular-winglet pair for identical parameters. 

e Maximum heat transfer for the case with rectangular winglet 
pair occurs at attack angle range between 45° and 65° and this 
maximum is smaller than that of delta-winglet pair of VG. 
While at 90° a sharp decrease in Nu number was observed. 

e Combined ribs and delta-winglet vortex generator give sig- 
nificant increase in heat transfer rate over the smooth channel. 

e A dramatic enhancement of heat transfer was observed when 
winglets placed on louvered of louvered fin heat exchanger. 

e The position of VG plays an important role on the performance 
of heat transfer and Nu number decreases with the space 
between the LVG pair decreases. 

e Longitudinal vortex generators can improve the synergy 
between the velocity and temperature field not only in the 
region near LVG but also in the large downstream region 
of LVG. 

e Heat transfer becomes steady at lower Re number while it 
becomes unsteady at the higher Re number values when VG 
is used. 

e The heat transfer rate and pressure drop directly increase with 
angle of attack, mass flow rate, length of vortex generator, and 
height of vortex generator (which should not be more than 
half of the channel height). 

e VGs were successfully used in narrow and microchannels to 
improve the heat transfer rate. 


This paper also presents an overview of recent developments 
in the experimental results and the numerical correlations on 
thermal conductivity using nanofluids. From literature, it can be 
summarized that nanofluids have high thermal conductivity 
compared to the conventional base fluids which is affected by 
several parameters. Nanoparticle concentration, nanoparticle size, 


nanoparticle type, nanoparticle shape, nanofluid temperature, 
base fluid type, and thermal stability are the parameters that 
influence the thermal conductivity of nanofluids. According to the 
high value of thermal conductivity of nanofluids, nanofluids are 
very suitable to be used as a working fluid in a wide range of 
applications especially when the pressure drop penalty accom- 
plished the use of nanofluid is too small or nil. As a conclusion: 


e A small amount of nanoparticles suspended in a base fluid 
provides substantially higher thermal conductivity than that of 
base fluid. 

e In order to find out the temperature dependence of thermal 
conductivity of nanofluids, experiments were performed at 
room temperature up to 420K. No paper has been yet 
published to report experiments dealing with low range of 
temperatures ( < 0 °C). 

e For forced convection heat transfer, convective heat transfer 
coefficient increases with nanoparticle concentration when 
the mass flow rate is kept constant. 

e Theoretically, viscosity, as well as the effective thermal con- 
ductivity of nanofluid, plays a key role in the natural convec- 
tion heat transfer. 

e Distinct deviations in the values of thermal conductivities of 
nanofluids obtained by models. Lack of theoretical under- 
standing of the heat transfer mechanism using nanofluid. 
Therefore, further investigations are too necessary in future. 
Lack of agreement in the experimental results from different 
groups about the pressure drop associated to the suspension of 
nanoparticles in the base fluid. 


Nanofluids open wide directions in the applications of thermal 
systems. Therefore, often research efforts are necessary to focus 
on the applications of nanofluids such as engineering, medical 
and space applications. 


6. Outlook and future challenges 


Further efforts are too necessary in future to focus on finding 
out new models and correlations to predict precise thermal 
conductivities with small deviation with the experimental results. 
Furthermore, future work should be done to give an explanation 
of heat transfer mechanism of nanofluids. The development in 
nanotechnology will be fruitful and helpful for nanofluid investi- 
gation. Combined VG and nanoparticles are new innovative idea 
in order to improve the performance of thermal systems which 
did not reported yet in the open literature. Wide opportunities 
will be opened when authors focus their efforts on using nano- 
fluids combined with VGs in micro-scale heat exchangers and 
heat sinks. The results of the heat transfer and pressure drops of 
nanofluids in different shapes of compact heat exchangers and 
heat sinks gave high deviations compared to the results obtained 
by the theoretical correlations and experimental results obtained 
by using base fluids. So the behavior of the nanofluids is definitely 
different than the trend of the base fluids which motivate the 
investigators to concentrate their researches to perform a lot of 
experiments to study the characteristics of nanofluids in several 
micro-scale channels configurations to derive correlations for 
Nusselt number and friction factor particular for nanofluids to 
aid authors to compare the numerical data with the experiments. 
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